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Preface 


This  paper  provides  a  summary  of  a  series  of  wind  wave  frequency- 
direction  spectral  observations  made  with  a  two-dimensional,  high-resolution 
directional  wave  gage  of  opportunity.  These  observations  were  the  result  of 
data  shared  by  researchers  from  Scripps  Institution  of  Oceanography  (SIO)  in 
execution  of  the  Sources  of  Ambient  Micro-Seismic  Noise  (SAMSON)  experiment. 
The  work  was  motivated  by  an  opportunity  to  conduct  a  360-deg  analysis  of  data 
from  a  24 -gage  array  and  thus  find  an  extraordinary  observation  set  of 
directionally  distributed  wave  energy  to  augment  understanding  and  modeling  of 
nearshore  processes  that  affect  coastal  engineering  projects.  This  effort  was 
authorized  by  Headquarters ,  US  Army  Corps  of  Engineers  (HQUSACE) ,  under  Civil 
Works  Research  Work  Unit  32484,  "Directionality  of  Waves  in  Shallow  Water," 
Coastal  Flooding  Program.  Funds  were  provided  through  the  Coastal  Engineering 
Research  Center  (CERC) ,  US  Army  Engineer  Waterways  Experiment  Station  (WES) , 
under  the  program  management  of  Ms.  Carolyn  M.  Holmes,  CERC.  Messrs.  John  H. 
Lockhart,  Jr.;  John  G.  Housley;  Barry  W.  Holliday;  and  David  A.  Roellig  were 
HQUSACE  Technical  Monitors. 

This  data  summary  was  prepared  by  Dr.  Charles  E.  Long  at  CERC’s  Field 
Research  Facility  (FRF)  in  Duck,  NC,  under  the  direct  supervision  of 
Mr.  William  A.  Birkemeier,  Chief,  FRF,  and  Mr.  Thomas  W,  Richardson,  Chief, 
Engineering  Development  Division,  CERC;  and  under  the  general  supervision  of 
Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Director  and  Assistant 
Director,  CERC,  respectively. 

The  SAMSON  experiment  was  planned  by  Drs,  Robert  T.  Guza,  SIO,  and  Joan 
M.  Oltman-Shay  of  Quest  Integrated,  Inc.,  Kent,  WA.  The  directional  array  was 
designed  by  Drs.  Guza  and  Thomas  H.  C,  Herbers,  SIO,  Data  processing  software 
was  written  by  Dr.  Oltman-Shay.  The  concept  of  continued  data  collection 
after  the  primary  SAMSON  experiment  is  credited  to  Dr.  C.  Linwood  Vincent, 
CERC.  The  efforts  of  Ms.  Janean  Shirley,  Information  Technology  Laboratory, 
in  editing  this  report  are  greatly  appreciated. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr. 

Robert  W.  Whalin.  Commander  was  COL  Leonard  G.  Hassell,  EN. 
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1.  During  October  and  early  November  1990,  an  experiment  involving  the 
detection  of  interacting  wave  trains  was  conducted  at  the  Coastal  Engineering 
Research  Center  Field  Research  Facility  (FRF).  Known  as  SAMSON  (Sources  of 
Ambient  Micro-Seismic  Ocean  Noise),  the  experiment  was  designed  to  detect 
relatively  unattenuated  near-bottom  pressure  signals  arising  from  nonlinear 
interactions  of  wind  wave  trains  traveling  in  opposing  directions.  These 
small  but  detectable  pressure  signals  are  believed  to  initiate  ocean  bottom 
seismic  waves  that  can  propagate  long  distances  and  appear  frequently  in 
seismograms  worldwide.  In  detecting  both  the  presence  of  interacting  wave 
trains  and  the  resulting  pressure  signals,  results  of  the  SAMSON  experiment 
can  be  used  to  test  hypotheses  concerning  the  initiation  of  such  seismic  waves 
near  a  coast  where  nearshore  bathymetry  can  act  as  a  partial  reflector  of 
ocean  wind  waves  and  thereby  induce  the  requisite  opposing  wave  trains.  The 
experiment  was  designed  and  executed  under  the  direction  of  Drs.  Robert  T. 

Guza  and  Thomas  H.  C.  Herbers  of  the  Center  for  Coastal  Studies  at  Scripps 
Institution  of  Oceanography,  University  of  California,  San  Diego,  and  Dr.  Joan 
M.  Oltman-Shay  of  Quest  Integrated,  Inc.,  Kent,  Washington,  with  the  experi¬ 
ment  site,  staging  area,  and  logistical  support  provided  by  the  FRF. 

2.  The  primary  instrument  in  the  SAMSON  experiment  was  a  two-dimension¬ 
al  array  of  24  bottom-mounted  pressure  gages  located  in  a  nominal  water  depth 
of  13  m.  By  design,  this  depth  is  deeper  than  is  normally  used  for  measure¬ 
ment  of  ocean  wind  waves.  The  idea  was  to  let  the  water  column  act  as  a 
filter  for  pressure  signals  arising  directly  from  wind  waves  above  a  certain 
cutoff  frequency.  In  that  way,  any  detectable  pressure  signal  above  the  cut¬ 
off  frequency  could  be  interpreted  as  being  from  the  wave-wave  interactions 
sought  in  the  experiment.  At  frequencies  below  the  cutoff,  the  pressure  gages 
could  readily  detect  wind  wave  signals.  As  a  result,  the  array  of  gages  was 
an  excellent  directional  detector  for  waves  in  the  lower  part  of  the  wind  wave 
frequency  band,  from  about  0.04  to  0.24  Hz. 
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3.  The  SAMSON  pressure  gage  array  was  deployed  in  September  -  90  and 
the  formal  experiment  took  place  in  October  and  early  November.  Scheduling 
recovery  of  the  gages  immediately  after  the  experiment  was  deemed  unwise 
because  of  the  onset  of  winter  and  the  high  likelihood  of  poor  working 
conditions.  Hence,  the  gages  were  left  in  place  until  the  following  spring 
and,  because  it  was  a  fully  functioning  instrument  s  t,  data  were  collected 
for  as  long  as  the  gages  continued  to  function,  though  this  was  not  expected 
to  be  for  very  long,  given  the  harshness  of  the  environment  in  which  they  were 
deployed.  Curiousl> ,  22  of  the  24  gages  continued  to  function  perfectly  until 
the  first  part  of  June  1991,  at  which  time  the  gages  were  recovered.  Data 
from  this  post-experiment  phase  of  instrument  deployment,  mid-November  1990  to 
the  end  of  May  1991,  have  been  made  available  to  the  FRF  for  directional 
analysis.  This  report  addresses  the  nature  and  results  of  that  analysis. 

4.  The  geometry  of  the  SAMSON  array  conformed  to  guidelines  given  both 
by  Davis  and  Regier  (1977)  for  a  method  of  wave  directional  detection  known  as 
Maximum  Likelihood  Estimation  (MLE)  and  for  an  improvement  on  this  method, 
known  as  Iterative  Maximum  Likelihood  Estimation  (IMLE),  derived  by  Pawka 
(1982,  1983).  Early  attempts  to  process  the  SAMSON  data  with  the  MLE  method 
were  less  than  satisfactory  for  some  aspects ,  primarily  the  estimation  of 
reflection  coefficients,  of  the  analysis.  Consequently,  the  IMLE  method  was 
employed  and  it  provided  much  better  results. 

5.  Because  the  SAMSON  array  is  two-dimensional,  it  can  detect  wave 
energy  propagating  in  any  direction.  This  allows  estimation  of  energy 
propagating  seaward  and  avoids  the  180-deg  ambiguity  present  in  analysis  of 
data  from  the  FRF  linear  array  directional  wave  gage  reported  previously  (Long 
and  Oltman-Shay  1991,  Long  1991a,  1991b).  In  this  report,  seaward-propagating 
energy  is  called  'reflected"  energy  because  reflection  of  incident  waves  from 
nearshore  bathymetry  is  its  most  likely  cause,  though  other  processes  may 
possibly  be  present.  Determination  of  reflected  directional  spectra  is  of 
enormous  interest  in  coastal  research  because  of  the  myriad  processes  that  are 
possible  theoretically  if  significant  reflection  is  present.  Hence,  this 
report  is  the  first  to  describe  a  well-resolved  estimation  of  the  full  360-deg 
wave  climate  at  the  FRF. 

6.  The  beginning  text  of  this  document  describes  and  clarifies  the 
substantial  information  contained  in  the  appendixes .  A  brief  overview  is 
given  of  the  measurement  site,  instrumentation,  data  collection,  and  method  of 


directional  spectral  estimation.  These  subjects  are  described  in  other 
publications  to  which  the  reader  is  referred  for  greater  detail.  Following 
the  overview  is  a  description  of  the  archived  frequency-direction  spectra  and 
some  characterizing  bulk  parameters  that  can  be  derived  from  them.  Appendix  A 
is  a  listing  of  these  characterizing  parameters  and  is  intended  to  be  used  as 
a  kind  of  catalog  of  the  set  of  spectra.  Appendix  B  contains  graphs  of  time 
series  of  some  of  these  parameters  as  a  pictorial  augmentation  of  the  informa¬ 
tion  in  Appendix  A.  Appendix  C  describes  tests  of  the  ability  of  the  SAMSON 
array  together  with  the  IMLE  algorithm  to  detect  low  levels  of  reflected  wrve 
energy  in  the  presence  of  higher  levels  of  incident  wave  energy  and  varying 
levels  of  background  noise.  It  provides  a  measure  of  confidence  for  the 
reflection  coefficients  reported  in  Appendixes  A  and  B. 

Measurement  Site 

7.  As  shown  in  Figure  1,  the  FRF  is  located  on  the  barrier  island  chain 
of  coastal  North  Carolina.  A  detailed  description  of  the  layout,  function, 
and  capabilities  of  the  FRF  is  given  by  Birkemeier  et  al.  (1985).  A  detailed, 
quantitative  description  of  the  climate  at  the  FRF,  as  determined  from  its 
arsenal  of  instrumentation  and  for  the  11  years  of  its  existence,  is  given  by 
Leffler  et  al.  (1992).  Of  particular  relevance  to  directional  wave  studies 
are  the  wave-steering  bathymetry  and  wave -generating  winds. 

8.  As  regards  the  former,  the  coastline  in  the  vicinity  of  the  FRF  is 
nearly  straight  for  several  tens  of  kilometers  north  and  south  (Figure  1).  It 
is  oriented  such  that  a  shore-normal  line  (directed  seaward)  is  very  nearly  70 
deg  from  true  north.  Waves  and  onshore  winds  can  approach  this  site  along  an 
easterly  180-deg  arc  from  340  deg  to  160  deg  true.  The  adjacent  continental 
shelf  is  wide,  relatively  shallow,  and  of  somewhat  complex  bathymetry.  The 
direction  of  nearest  approach  of  the  100-m  isobath,  which  indicates  the  shelf 
break,  is  10  deg  to  15  deg  south  of  east  and  is  about  80  km  distant.  A 
typical  bottom  slope  for  the  shelf  is  1  m/km,  but  this  is  interrupted  by 
numerous  features  of  1-  to  10-km  horizontal  scales  and  10-m  vertical  scales 
scattered  irregularly  across  the  shelf. 

9.  Within  a  few  kilometers  of  the  FRF,  the  offshore  bathymetry  is  more 
regular,  with  isobaths  nearly  shore-parallel  and  a  bottom  slope  of  about 

2  m/km  (Figure  2).  Some  irregularities  exist.  Within  about  300  m  of  the 


* 


*1 


•  • 


6 


7 


Figure  2.  Approximate  nearshore  bathymetry  and  SAMSON  array 
location  (seawardmost  set  of  dots) 

shore,  there  exists  a  complex  and  mobile  bar  system  (Birkemeier  1984). 
Furthermore ,  waves  and  currents  have  created  some  irregular  bathymetry  in  the 
vicinity  of  the  FRF  research  pier,  which  extends  about  600  m  offshore  (Miller, 
Birkemeier,  and  DeWall  1983). 

10.  The  site  is  subject  to  a  variety  of  climates  that  give  rise  to  a 
diverse  set  of  directional  wave  conditions.  Primary  sources  of  high-energy 
waves  are  winds  associated  with  hurricanes  and  frontal  passages.  Though 
Hurricane  Lily  moved  up  the  middle  Atlantic  Ocean  during  the  formal  SAMSON 
experiment,  no  other  hurricanes  affected  the  FRF  wave  climate  in  the  period 
covered  by  this  report.  Low-pressure  weather  fronts,  of  which  several  crossed 
the  FRF  site  during  this  reporting  period,  were  typically  oriented  northeast- 
southwest,  with  strong  wave -generating  winds  coming  from  the  northeast. 


8 


11.  As  a  consequence  of  the  complex  offshore  bathymetry  and  the  close 
proximity  of  wave -generating  wind  systems,  it  is  expected  that  directional 
wave  spectra  will  exhibit  a  rather  broad,  complex  structure.  This  expectation 
has  been  realized  in  studies  at  the  FRF  using  a  linear  array  directional  wave 
gage  deployed  about  900  m  offshore  in  8  m  of  water  (Long  and  Oltman-Shay 
1991).  With  the  exception  of  well-sorted  waves  from  distant  storms  and  waves 
in  the  initial  growth  stages  under  strong  steady  local  winds,  most  spectra 
were  found  to  be  directionally  broad.  Because  of  the  180-deg  ambiguity 
inherent  in  a  linear  array,  it  is  possible  that  part  of  this  directional 
spread  could  be  due  to  reflected  wave  energy,  which  is  folded  over  to  appear 
as  incident  wave  energy  in  the  linear  array  data  processing  algorithm.  The 
ability  to  perform  a  360-deg  wave  analysis  with  the  SAMSON  array  helps  to 
ascertain  the  extent  of  the  error  induced  by  the  180-deg  ambiguity  of  the 
linear  array,  in  addition  to  estimating  the  amount  of  seaward- propagating 
energy . 


Instrumentation 


12.  The  primary  directional  wave  detector  used  in  this  study  was  a  q 

two-dimensional  array  of  wave  gages  coupled  with  the  high- resolution  IMLE 

algorithm  fo-_  conversion  of  gage  data  to  directional  spectra.  The  array  of 
gages  sampled  the  near-bottom  pressure  field  and,  through  linear  wave  theory, 
provided  estimates  of  sea-surface  displacement  at  several  points  in  (hori-  ^ 

zontal)  space.  The  mathematical  treatment  of  these  data  to  obtain  estimates 
of  wave  directionality  involves  two  tuning  procedures  to  optimize  results. 

One  procedure  is  the  use  of  subsets  of  gages  to  reduce  the  effective  size  of 

the  array  for  waves  with  short  wavelengths.  The  other  procedure  is  the  • 

elimination  of  select  gages  from  analysis  in  the  event  that  the  array  provides 

too  much  redundant  information  in  a  particular  observation. 

BaSlS-, ftEIfllE 

13.  The  SAMSON  array  consisted  of  24  pressure  gages  buried  about  0.1  m  • 

below  the  bottom  in  the  vicinity  of  the  13-m  isobath  about  1,700  m  offshore 

and  slightly  to  the  north  of  the  research  pier.  This  is  shown  in  Figure  2  as 
the  seawardmost  set  of  dots.  In  addition  to  the  design  requirements  mentioned 
above,  its  location  satisfies  three  constraints.  First,  the  array  is  general-  • 

ly  outside  the  surf  zone  so  that  linear  wave  theory  is  applicable  in  data 
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processing.  Second,  it  is  in  water  shallow  enough  that  signals  from  4-sec 
waves  are  detectable  above  background  noise  at  the  bottom-mounted  gages. 
Third,  the  array  is  located  away  from  the  irregular  isobaths  around  the  pier 
and  in  the  nearshore  bar  system,  which  helps  minimize  bathymetrically  induced 
inhomogeneities  in  the  wave  field. 

14.  Figure  3  is  an  enlarged  view  of  the  array  layout  and  shows  gage 
spacing  as  well  as  the  gage  numbering  scheme.  Spacing  between  the  gages  in 
the  array  corresponds  reasonably  well  to  the  array-design  criterion  posed  by 
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Davis  and  Regier  (1977)  that  every  gage  pair  have  a  unique  separation  or 
spatial  lag.  There  are  a  large  number  of  redundant  lags  built  into  the  array, 
so  that  multiple  estimates  of  cross -spectral  densities  can  be  made  to  check 
for  spatial  inhomogeneity  in  the  wave  field  and  to  reduce  system  noise.  There 
are  more  unique  lags  with  which  to  fill  in  the  spatial  cross -spectral  matrix 
from  which  wave  direction  estimates  are  made.  Minimum  gage  spacing  is  4  m, 
maximum  spacing  (the  greatest  dimension  of  the  array)  is  256  m,  and  intermedi¬ 
ate  gage  spacings  are  in  multiples  of  4  m.  Unfortunately,  gages  12  and  14 
were  inoperative,  and  the  loss  of  gage  12  reduced  the  longshore  dimension  of 
the  array  to  192  m,  at  some  cost  to  direction-resolving  ability  for  the 
longest  (0.05-Hz)  waves  measured.  It  is  noted  that  gage  19  appears  to  be 
offline  from  the  basic  gage  pattern,  and  it  was,  in  fact,  slightly  misplaced. 
However,  this  does  not  make  it  ineffective.  In  IMLE  processing,  it  is  only 
necessary  for  the  true  gage  positions  to  be  known,  and  not  that  they  conform 
to  any  particular  pattern.  The  SAMSON  gage  positions  were  surveyed  to  an 
accuracy  of  better  than  1  m  after  the  gages  were  in  place.  At  worst,  this 
uncertainty  would  result  in  a  directional  error  of  about  1.8  deg  for  the 
shortest  (0.24-Hz)  waves  considered  and  about  0.6  deg  for  the  longest  (0.04- 
Hz)  waves. 

Sub -arrays 

15.  At  sufficiently  large  lags,  signals  from  a  natural  wave  field  tend 
to  become  uncorrelated.  As  a  result,  cross -spectra  between  widely  spaced 
gages  tend  to  become  mostly  noise.  At  very  small  lags,  correlations  and 
cross -spectra  are  difficult  to  discern  from  those  at  zero  lag  and  so  tend  to 
become  somewhat  redundant.  The  terms  large  and  small  are  generally  related  to 
the  wavelengths  A*  of  waves  under  consideration.  To  avoid  both  redundancy 
at  small  lags  and  the  problem  of  noise  at  large  lags,  the  SAMSON  array  data 
were  processed  by  using  subsets  of  gages,  or  sub-arrays,  based  on  ranges  of 
wavelengths  being  considered. 

16.  Figure  4  illustrates  the  sub-arrays  and  the  ranges  of  wave  frequen¬ 
cy  f  used  with  each  sub- array  in  the  present  analysis.  Note  that  gages  12 
and  14,  from  which  reliable  data  were  not  obtained,  are  not  shown.  For  the 
longest  waves,  all  the  working  gages  were  used  except  those  that  introduced 


S> 

*) 


•  •! 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  D) . 


11 


Figure  4.  Sub-arrays  and  their  associated  frequency  bands 

the  shortest  (4-a)  lags,  i.e.,  gages  5  and  22.  This  left  20  gages  for 
processing  data  in  the  frequency  range  0.05  <  f  (Hz)  <0.10  or  wavelength 
range  103  <  X  (a)  <  220  (assuming  a  nominal  water  depth  of  13  a) ,  and 
included  the  longest  possible  lags. 

17.  For  shorter  wavelengths,  the  selection  criterion  was  that  the 
longest  dimension  of  the  sub -array  not  exceed  about  three  wavelengths  of  the 
shortest  waves  in  a  given  band.  For  the  low  aid- range  of  frequencies, 

0.10  <  f  (Hz)  <0.16  ,  corresponding  to  wavelengths  55  <  X  (a)  <  103  ,  this 
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procedure  eliminated  gage  19  from  analysis,  leaving  19  gages  in  the  sub-array. 
For  the  high  mid-range  of  frequencies,  0.16  <  f  (Hz)  <0.20  ,  or 
38  <  X  (m)  <  55  ,  gages  1,  2,  13,  and  21  were  eliminated,  but  the  short-lag 
gages  5  and  24  were  included,  for  a  total  of  17  gages.  For  the  highest 
frequencies.  0.20  <  f  (Hz)  <  0.24  ,  27  <  \  (m)  <  38  ,  gages  3,  11.  and  20 
were  excluded  from  the  previous  set,  leaving  14  gages. 

aa.ta.:^flE£l,yg....gag6.  ,.g  1  imiaat Ian 

18.  As  noted  by  Davis  and  Regier  (1977),  wave  fields  may  occur  for 
which  the  cross -spectral  matrix  of  a  given  array  geometry  contains  nearly  re¬ 
dundant  information.  When  this  happens,  it  becomes  impossible  to  compute  the 
inverse  of  the  cross -spectral  matrix,  a  process  that  is  central  to  the  IMLE 
algorithm,  and  no  directional  information  can  be  found.  This  condition  can 
occur  when  the  wave  field  has  a  particularly  simple  directional  distribution, 
such  as  a  nearly  unidirectional  wave  train,  or  when  the  wave  field  is  more 
complex,  but  crosses  the  array  of  gages  in  such  a  way  that  many  pairs  of  gages 
have  approximately  equal  cross -spectral  values.  This  condition  is  more  likely 
to  occur  when  the  number  of  gages  is  large,  as  in  the  present  case.  One 
solution  to  this  dilemma  is  to  redeploy  the  gages  in  a  pattern  that  does  not 
yield  so  much  redundant  information,  but  such  a  procedure  is  not  likely  to  be 
successful  at  a  site,  like  the  FRF,  where  the  wave  field  is  known  to  have 
quite  diverse  characteristics.  In  any  case,  such  a  procedure  was  not  practi¬ 
cal  in  the  present  study  because  gage  positions  were  fixed  before  the  study 
began. 

19.  The  approach  used  here  was  to  eliminate  the  data  of  one  or  more 
gages  from  analysis  and  try  the  directional  computation  again  until  success 
was  achieved.  The  order  of  gage  elimination  was  assigned  a  priori  by  examin¬ 
ing  the  nature  of  coverage  provided  by  a  reduced  array  when  transformed  into 
lag  coordinates.  For  the  sub-arrays  associated  with  long  waves,  gages  were 
first  eliminated  that  gave  small  lags  or  small  lag  increments,  but  still 
evenly  covered  the  largest  possible  spans  in  both  dimensions.  A  second 
criterion,  to  avoid  aliasing,  was  that  any  remaining  gages  not  have  lag 
Increments  larger  than  half  of  the  longest  wavelength  to  be  analyzed.  For 
short-wave  sub-arrays,  gages  were  first  eliminated  that  gave  large  lags.  The 
second  criterion  for  short  waves  was  the  same  as  for  long  waves.  Table  1 
lists  the  gages  used  in  each  frequency  range  and  the  numerical  order  in  which 
they  were  to  be  eliminated.  For  example,  for  the  lowest  frequencies,  gage  24 
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Table  1 


Order  of  Gage  Elimination  from  SAMSON  Sub-Arravs 


Frequency  range:  0.04  <  f  (Hz)  <0.10 

Sub -array  gages: 

1.  2,  3,  4.  6.  7.  a.  9.  10.  11.  13,  15,  16.  17,  10,  19,  20.  21.  23.  24 

Elimination  order: 

24,  6.  16,  a.  23,  10.  19.  11,  1,  13,  21,  2.  20,  IB,  4,  9,  17.  15,  3,  7 


Frequency  range: 

0.10 

< 

f  (Hz)  <  0.16 

Sub -array  gages: 

1,  2,  3. 

4,  6,  7,  8. 

9, 

10,  11.  13,  15.  16, 

17, 

18. 

20. 

21, 

23.  24 

Elimination  order 

: 

24,  1,  2, 

13,  4,  8,  23,  1 

LI. 

21,  20,  3,  15,  10, 

18, 

7. 

16, 

17, 

6,  9 

Frequency  range: 

0.16 

< 

f  (Hz)  <  0.20 

Sub -array  gages: 

3.  4, 

3.  8,  7.  8, 

9. 

10.  11,  15.  16,  17, 

18. 

20, 

22. 

23. 

24 

Elimination  order 

: 

20,  11, 

4,  10,  24.  8. 

3. 

13,  18,  16,  22,  5. 

6, 

23. 

17, 

7, 

9 

Frequency  range:  0.20  <  f  (Hz)  <0.24 

Sub -array  gages: 

4.  5.  6.  7,  a,  9,  10,  15.  16,  17,  18,  22,  23,  24 

Elimination  order: 

18,  10,  17,  08,  15.  16,  9,  24,  7,  22,  5,  23.  4,  6 


Note:  Gages  12  and  14  are  excluded. 


is  eliminated  first,  followed  by  gage  6,  gage  16,  gage  8,  etc.  Though  a 
complete  elimination  order  is  listed,  a  limit  was  set  that  at  least  four  gages 
must  remain  for  a  directional  estimate  to  be  made.  This  limit  was  not  met  in 
practice,  though  the  system  was  occasionally  reduced  to  five  or  six  gages. 
Pressure  gazes  and  data  collection  hardware 

20.  Each  pressure  gage  is  a  Setra  Systems  Model  280E  capacitive -type 
transducer  having  a  range  of  0  to  50  psia*.  The  manufacturer's  specifica¬ 
tions  indicate  an  accuracy  of  0.11  percent  of  full  scale,  or  the  pressure 


*  A  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI  (metric) 
units  is  presented  on  page  3. 
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equivalent  of  ±0.037  a  of  water.  Analog  pressure  signals  were  collected  in  an 
underwater  terminal  where  they  were  discretized  to  12-bit  binary  form, 
multiplexed,  and  encoded  for  transmission  down  a  single  data  line.  Digitiza¬ 
tion  of  the  signal  to  12-bit  binary  form  resulted  in  a  discretization  uncer¬ 
tainty  of  the  equivalent  of  0.002  m  of  water.  At  the  landward  end  of  the  data 
line,  the  encoded  signal  was  decoded  and  sent  to  a  parallel  port  on  a  personal 
computer  for  data  acquisition.  During  the  post-experiment  period,  data  files 
from  the  personal  computer  were  periodically  transferred  to  the  higher 
capacity  hard  disks  of  the  FRF’s  Digital  Equipment  Corporation  VAX  11/750 
computer  for  accumulation  and  eventual  archiving  on  magnetic  media. 

Data  Collection 

21.  Signals  from  each  of  the  24  pressure  gages  were  sampled  at  4  Hz  and 
stored  in  records  of  768  samples.  A  total  of  1,280  records,  or  10,240  data 
points,  were  collected  for  each  gage  for  an  elapsed  time  of  2  hr  50  min  40  sec 
for  each  collection.  Starting  times  for  collections  during  the  main  SAMSON 
experiment  were  daily  at  0100,  0400,  0700,  1000,  1300,  1600,  1900,  and  2200 
Eastern  Standard  Time  (EST)  and  this  pattern  was  continued  into  the  post¬ 
experiment  period  until  1  December  1990.  At  that  time,  collection  start  times 
were  set  to  coincide  with  routine  FRF  observations  (Birkemeier  et  al.  1985), 
which  occur  daily  at  0100,  0700,  1300,  and  1900  EST.  This  collection  pattern 
was  maintained  until  the  gages  were  disconnected  for  retrieval  on  2  Jun  1991. 

22.  In  the  analysis  presented  here,  only  the  first  2  hr  16  min  32  sec 
(eight  records  of  4,096  points)  of  data  from  each  collection  was  used.  This 
procedure  resulted  in  data  sets  that  coincided  both  in  start  times  and 
duration  with  other  permanent  directional  wave  gages  at  the  FRF.  It  also 
allowed  consistency  in  processing  some  collections  that  were  terminated  early 
to  allow  time  for  file  transfers  from  the  data  collection  computer  to  the  FRF 
VAX  11/750  hard  disks. 


Data  Processing 

23.  Conversion  of  measured  time  series  to  estimates  of  frequency- 
direction  spectra  requires  products  of  frequency  spectral  estimates  from  all 
the  various  gages  in  the  array.  For  final  results  to  be  accurate,  raw  input 
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data  must  be  of  exceptionally  high  quality  so  that  spiky  or  drifty  data  from 
one  gage  do  not  contaminate  products  of  results  from  the  other  gages.  Hence, 
the  procedure  for  data  processing  is  to  check  the  raw  data  for  errors, 
estimate  the  frequency-direction  spectrum,  and  then  compute  some  bulk  para¬ 
meters  with  which  to  characterize  results . 

Error  checking 

24.  Because  multiple  gages  were  deployed  in  an  (assumed)  uniform  sea, 
certain  properties  of  raw  data  from  the  pressure  gages  should  be  identical.  A 
convenient  method  of  error  checking  is  to  compute  frequency  spectra  of  the  raw 
time  series  from  each  of  the  22  functioning  gages  and  overplot  the  results  on 
a  single  graph.  This  method  was  used  by  Drs.  Guza  and  Herbers  during  the  main 
experiment  and  has  also  been  incorporated  for  the  present  analysis.  Because 
of  the  depth  at  which  the  SAMSON  gages  were  deployed,  very  little  direct  wind 
wave  signal  is  expected  in  the  frequency  range  from  about  0.3  Hz  out  to  the 
Nyquist  frequency  of  2.0  Hz.  The  frequency  spectrum  in  this  frequency  band 
should  primarily  indicate  system  noise,  which  should  be  about  the  same  for 
each  gage.  Excessively  spiky  data  from  one  gage  will  raise  the  noise  level 
for  that  gage  relative  to  the  noise  levels  from  the  other  gages,  so  that  its 
signature  will  be  distinct  in  a  collection  of  overplotted  spectra.  Strong 
drifts  in  the  signal  can  appear  in  a  variety  of  ways  in  corresponding  spectra, 
but  will  generally  give  a  distinct  signal  at  low  frequencies.  In  the  pass 
band  of  wind  wave  frequencies  for  which  directional  estimates  are  computed 
(here  0.04  to  0.24  Hz),  one  expects  the  frequency  spectra  to  be  nearly 
identical . 

25.  Figure  5  is  an  example  of  one  set  of  overplotted  frequency  spectra. 
Semi -logarithmic  coordinates  have  been  used  and  pressure  gage  voltages  have 
been  converted  to  equivalent  heights  of  a  static  water  column,  but  no  other 
transformations  have  been  used.  As  can  be  seen,  signals  in  the  wind  wave 
frequency  pass  band  are  all  very  nearly  alike,  indicating  that  all  gages  are 
functioning  reasonably  well.  The  noise  floor  at  high  frequencies  is  very  low 
relative  to  the  wind  wave  signal  and  is  nearly  uniform  for  all  gages.  The 
spikes  in  the  spectra  at  high  frequencies  are  due  to  low  levels  of  commercial 
60-Hz  power  line  noise  present  in  the  system.  The  inset  graph  in  Figure  5 
identifies  the  gages  used  and  shows  the  deviations  of  their  mean  values  from 
the  median  of  all  the  mean  values  after  adjusting  for  their  relative  depths. 
Also  shown  is  the  deviation  of  mean  water  level  from  still -water  level  (based 
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SAMSON  Array  Frequency  Spectra  (Bottom) 
Date:  05  Dec  90  Time:  0700 

Missing  gages:  NONE 


Figure  5.  Example  of  overplocted  frequency  spectra 


on  the  1929  National  Geodetic  Vertical  Datum)  so  that  tidal  level  and  tidal 
variation  can  be  ascertained.  If  a  particular  gage  becomes  more  noisy  than 
the  others,  it  can  then  be  readily  identified. 

26.  Tests  with  this  kind  of  checking  routine  indicated  that  the  SAMSON 
gages  and  data  collection  system  were  extremely  reliable.  At  no  time  during 
the  post -experiment  period  were  any  of  the  22  functioning  gages  found  to  have 
any  significant  errors. 

Frequency-direction  spectra 

27.  Estimation  of  frequency-direction  spectra  is  done  in  four  parts. 
First,  time  series  of  pressure  data  from  each  gage  are  Fourier  transformed  to 
the  frequency  domain.  Second,  these  transforms  are  converted  to  sea- surface 
displacement  transforms.  Third,  cross  spectra  of  sea-surface  displacement  are 
computed  for  all  unique  gage  pairs  at  each  frequency.  Finally,  an  estimate  is 
made  of  a  directional  distribution  of  wave  energy  that  corresponds  to  the 
computed  spatial  variation  in  cross -spectral  density  for  each  frequency. 

28.  The  Fourier  transform  is  conventional.  A  32, 768 -point  time  series 
is  divided  into  15  half -overlapping  segments  of  4,096  points.  Segments  are 
tapered  with  a  Kaiser-Vessel  window  (a  modified  Bessel  function  of  the  first 
kind,  compensated  uniformly  for  loss  of  variance  due  to  windowing)  and  fast 
Fourier  transformed.  An  intermediate -resolution  transform  is  found  by 
averaging  the  15  transformed  segments,  frequency  by  frequency.  Final  trans¬ 
forms  are  found  by  then  averaging  results  over  10  adjacent  frequency  bands. 
Final  resolution  bandwidth  is  0.00976  Hz,  and  degrees  of  freedom  are  at  least 
150  (assuming  eight  contiguous  segments  and  ignoring  any  gain  from  lapped  seg¬ 
ments).  Transform  estimates  are  retained  for  20  frequency  bands  with  band- 
center  frequency  ranging  from  0.054  to  0.240  Hz. 

29.  Conversion  of  pressure  signals  at  depth  to  water-surface  displace¬ 
ment  is  done  through  the  linear  wave  theory  pressure  response  factor  as 
described  in  the  Shore  Protection  Manual  (SPM)  (1984).  After  this  conversion, 
complex  cross  spectra  in  the  form  of  coincident  and  quadrature  spectra  are 
computed  in  the  conventional  way  (Bendat  and  Piersol  1971,  Jenkins  and  Watts 
1968)  for  all  unique  gage  pairs.  Cross -spectral  estimates  at  a  given  frequen¬ 
cy  are  then  stored  in  a  matrix  of  complex  numbers,  the  indices  of  which  can  be 
used  to  address  data  arrays  containing  gage  position  for  the  computation  of 
gage  separation  distance,  or  lag  space,  in  preparation  for  directional 
spectral  estimation  at  that  frequency. 
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30.  As  mentioned  in  the  introduction,  conversion  of  cross -spectral 
patterns  in  lag  space  to  directional  spectra  is  done  with  the  IMLE  algorithm 
described  by  Pawka  (1983).  The  algorithm  is  also  described  in  application  to 
data  from  heave-pitch-roll  buoys  by  Oltman-Shay  and  Guza  (1984).  Specifical¬ 
ly,  data  were  processed  using  Equations  1  to  4  of  Pawka 's  (1983)  paper  with  a 
convergence  exponent  of  1.0,  a  convergence  coefficient  of  5.0,  and  a  maximum 
of  30  iterations. 

31.  The  algorithm  used  here  made  a  spectral  estimate  for  every  1  deg  of 
azimuthal  direction.  These  were  smoothed  over  2 -deg  arcs  to  create  180  bins 
covering  a  complete  circle  with  which  to  describe  the  directional  distribution 
of  wave  energy. 

32.  The  primary  result  of  data  processing  is  an  estimate  of  the 
discrete  frequency-direction  spectrum  S(fn,0B)  ,  which  represents  the  vari¬ 
ance  of  sea- surface  displacement  per  frequency  resolution  bandwidth  df 

(-  0.00976  Hz)  per  direction  resolution  arc  d 6  (-  2  deg),  where  fn  is  the 

nth  of  N  -  21  discrete  frequencies  and  8m  is  the  mch  of  M  -  180 
discrete  directions.  In  this  work,  direction  is  considered  to  be  the  angle 
from  which  wave  energy  is  coming,  measured  counterclockwise  from  shore-normal. 
It  is  noted  that  parametric  definitions  below  are  given  in  terms  of  these 
discrete  definitions  rather  than  in  terms  of  the  continuous  functions  they 
approximate . 

33.  Numerical  values  of  S(fn,0B)  can  range  over  many  orders  of 
magnitude,  depending  on  the  amount  of  energy  in  a  given  frequency  band  and 
direction  arc,  and  this  can  require  space -consuming  formats  for  archiving 
data.  To  simplify  this  problem,  frequency-direction  spectra  can  be  saved  in 
the  form  of  directional  distribution  functions  D(fn,0o)  defined  by 


s(f„) 


(1) 


where  S(fn)  is  the  frequency  spectral  density  at  frequency  fn  .  The 
directional  distribution  function  has  units  of  deg"1,  and  its  integral  with 
respect  to  direction  over  all  directions  is  unity. 

34.  The  frequency  spectrum  in  Equation  1  represents  the  total  over  all 
directions  of  sea- surface  variance  per  frequency  bandwidth  and  is  defined  in 
terms  of  the  frequency-direction  spectrum  by 


•  • 
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(2) 


M 


s<fn>  -  £  s<f„,*.)  at 

B»1 


where  M  -  180  as  mentioned  above.  Note  that  this  is  identical  to  a  conven¬ 
tional  frequency  spectrum  that  would  result  from  a  time  series  of  sea- surface 
displacements  at  a  single  point  in  space.  Because  it  is  an  integral  of  the 
frequency-direction  spectrum,  it  is  called  the  integrated  frequency  spectrum. 

35.  A  directional  analog  of  the  frequency  spectrum  is  the  integrated 
direction  spectrum,  found  by  summing  the  frequency-direction  spectrum  over  all 
frequencies  for  a  fixed  direction  arc.  It  is  computed  from 

K 

S(tf.)  «  £  S(fn.*.)df  (3> 

n- 1 

where  N  —  21  in  the  present  context.  Figure  6  shows  one  way  to  display  the 
frequency- direction  spectrum  and  the  corresponding  integrated  frequency  and 
integrated  direction  spectra. 

Bulk  parameters 

36.  Several  parameters  have  been  computed  to  characterize  the  observed 
spectra.  Five  elementary  types  of  parameter  have  been  used  here:  water 
depth,  characteristic  wave  height,  peak  period  (or  its  inverse,  peak  frequen¬ 
cy),  peak  direction,  and  reflection  coefficient.  There  is  more  than  one  way 
to  define  some  of  these  parameters,  so  definitions  are  presented  here. 

37.  Water  depth.  Water  depth  is  found  from  the  median  of  the  mean 
depths  from  all  22  gages  in  the  array  after  adjusting  for  their  relative 
deployment  depths  and  applied  at  the  location  of  gage  1,  which  is  used  as  a 
reference  gage.  No  compensation  for  variations  in  atmospheric  pressure  was 
done.  Water  depth  represents  a  water  level  that  varies  with  low  frequency 
processes  like  tides.  Because  variations  in  mean  water  level  adjust  intersec¬ 
tions  of  the  water  surface  with  nearshore  bathymetry,  it  may  be  an  important 
variable  in  accounting  for  the  amount  of  observed  reflected  energy. 

38.  Characteristic  wave  height.  Characteristic  wave  heights  from 
spectral  observations  are  most  frequently  given  as  H^,  ,  which  is  four  times 
the  standard  deviation  of  sea -surface  displacement.  It  can  be  determined  from 
the  volume  under  the  frequency -direction  spectrum  by  the  equation 
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Figure  6 . 


Sample  displays  of  frequency  spectrum,  direction  spectrum, 
and  frequency-direction  spectrum 
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(4) 


N  M 

E  S(fn,em)dtdf 

n»l  ®«1 

It  can  be  found  equivalently  from  the  integrated  frequency  spectrum  by 


«L  -  16  £  S(fn)df 

n«X 


(5) 


which  is  its  conventional  definition,  or  from  the  integrated  direction 
spectrum  by 

hL  ■  16  E  S(6n)dt)  (6) 

m-1 

39.  Because  the  analysis  yields  a  360-deg  distribution  of  wave  energy, 
it  is  reasonable  to  distinguish  between  incident  wave  energy  and  outgoing  wave 
energy,  which  here  is  arbitrarily  called  reflected  wave  energy.  Characteris¬ 
tic  wave  heights  can  be  assigned  for  both  incident  and  reflected  energies,  and 
the  ratio  of  these  wave  heights  is  an  estimate  of  the  overall  reflection 
coefficient.  Such  a  ratio  is  rather  sensitive  to  noise  in  the  directional 
distribution  estimate  on  which  it  is  based,  so  an  attempt  is  made  here  to 
reduce  the  noise  level  for  parameters  involved  with  wave  reflection.  It  is 
noted  that  the  tails  of  the  integrated  direction  spectrum  in  Figure  6  do  not 
go  to  zero,  but  appear  to  level  off  at  «ome  finite  value  that  persists  over 
all  directions.  While  it  is  possible  that  some  wave  energy  propagates  in  all 
directions,  it  is  not  likely.  It  is  assumed,  therefore,  that  the  small 
persistent  level  of  the  spectrum  represents  the  accumulation  of  instrument, 
system,  and  analysis  noise  S„  .  In  this  report,  its  value  is  estimated  as 
the  average  of  the  10  smallest  values  of  the  discrete  integrated  direction 
spectrum  and  it  is  assumed  to  be  constant  for  all  directions. 

40.  A  characteristic  incident  wave  height  t  can  be  identified  from 
the  volume  under  the  frequency-direction  spectrum  for  waves  approaching  from 
directions  in  the  range  -TO  deg  <  8  <  90  deg  or,  equivalently,  the  area  under 
the  integrated  direction  spectrum  over  the  same  range  of  directions  after  the 
removal  of  noise.  Discrete  directions  in  the  present  analysis  are  defined  by 

6m  ■  180  deg  -  (a  -  1)  2  deg  (7) 

If  half  the  energy  in  the  90-deg  and  -90-deg  bins  and  all  the  energy  in  the 
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intervening  directional  bins  is  attributed  to  incident  energy  plus  noise,  a 
characteristic  incident  wave  height  is  the  sum  of  the  spectral  estimates 
reduced  by  the  estimated  noise,  or 

135 

hL,i  -  8[S<*46)  -  S„]d*  ♦  16  £  [S(«.)  -  SN]dtf  ♦  8  [S(ff1J6)  -  S„}d*  <»> 

n»«47 

In  a  similar  fashion,  if  the  other  half  of  the  energy  in  the  90-deg  and 
-90-deg  bins  and  the  remainder  of  the  offshore  propagating  energy  is  attribut¬ 
ed  to  reflected  energy,  a  characteristic  reflected  wave  height  H^,  r  is  given 

by 

Cr  =  16  E  IS(«.)  -  SJdtf  ♦  16[iS(*t#)  *  \S(8136)  -  sjdf 

m-1  l z  z 

(9) 

180 

*  16  £  [S(*B>  -  S„]d* 

b-137 

41.  Peak  frequency.  Peak  frequency,  which  has  the  generic  notation 
fp  ,  can  be  defined  at  least  two  ways.  One  way  is  to  find  the  frequency  (and 
direction)  at  which  the  frequency-direction  spectrum  is  maximum.  Another  way 
is  to  find  the  frequency  at  which  the  integrated  frequency  spectrum  is 
maximum.  This  is  the  more  conventional  definition  and  is  the  basic  definition 
used  here.  As  with  characteristic  wave  height,  a  distinction  is  made  between 
peak  frequencies  associated  with  incident  and  reflected  wave  energy.  The  peak 
frequency  of  the  incident  wave  spectrum  f  L  is  defined  as  the  frequency  of 
the  maximum  of  the  integrated  incident  frequency  spectrum  Sl(fn)  defined 
following  the  notation  above  for  H,*,  (and  without  removing  any  noise 
estimate)  as 

133 

W  *  ±s(fn,ett)de  ♦  £  s(fn,*jd«  *  is(fn,*1J6)d*  do) 


•  • 


In  a  similar  manner,  the  integrated  reflected  wave  spectrum  Sr(fn)  is  given 
by 

*5  180 

Sr(V  -  E  S(fn,*„)d 9  ♦  i[S(fB>*„)  ♦  s<fn,s136)]dtf  ♦  £  S(fn.#.)d*  (ID 
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However,  this  spectrum  was  not  used  to  find  the  reflected  peak  frequency 
because  of  the  propensity  for  the  tails  of  the  incident  peaks  to  cross  the 
shore -parallel  azimuths  (±90  deg)  and  dominate  the  usually  weaker  reflected 
energy.  For  example,  the  spectral  surface  in  Figure  6  crosses  the  90- deg 
azimuth  at  f  -  0.24  Hz  and  clearly  has  more  energy  than  anywhere  else  along 
the  reflection  arc,  even  though  that  surface  is  associated  with  the  incident 
spectrum.  To  avoid  this  effect  in  the  definition  of  the  reflected  peak 
frequency  fpT  ,  an  auxiliary  frequency  spectrum  S,(fn)  was  computed  using 
only  those  parts  of  S(f „,8m)  for  which  |®|  £  120  deg  ,  or 

30  ISO 

S.(fn)  -  £  S(fn,*B)dfl  ♦  £  S(fn,«B)d*  (12) 

m-1  m«152 

The  frequency  at  which  S,(fn)  is  a  maximum  is  then  identified  as  f p  r  . 

42.  The  distinction  between  fp>1  and  f r  may  be  important.  Unless 
there  are  very  strong  reflections,  the  incident  peak  frequency  is  virtually 
the  same  as  the  peak  frequency  one  would  obtain  from  a  time  series  at  a  single 
point  in  space.  If  the  beach  is  uniformly  reflective  (the  same  for  all 
frequencies),  reflected  peak  frequency  will  be  the  same  as  incident  peak 
frequency.  If  the  two  frequencies  are  not  equal,  the  beach  may  be  signifi¬ 
cantly  more  reflective  for  certain  frequencies  other  than  the  peak,  a  situa¬ 
tion  of  no  small  dynamic  importance  in  nearshore  processes. 

43.  Peak  period.  Peak  period  Tp  is  the  characteristic  wave  period 
associated  with  spectral  peak  frequency.  As  simply  inverses  of  peak  frequen¬ 
cies,  there  are  two  peak  periods  associated  with  a  directional  spectrum 
corresponding  to  the  two  peak  frequencies  defined  above.  The  incident  peak 
period  is  Tp>l  -  l/fp  t  and  the  reflected  peak  period  is  Tp  r  -  l/fp  r  . 

44.  Peak  direction.  Peak  direction  is  the  direction  at  which  the  most 
energy  occurs.  Given  the  generic  symbol  8p  ,  it  is  defined  as  the  peak  of 
the  direct: cnal  distribution  associated  with  a  peak  frequency.  As  with  the 
other  parameters  there  is  an  incident  peak  direction  8p  t  and  a  reflected 
peatc  direction  8pz  .  Incident  peak  direction  is  the  direction  of  the 
maximum  of  S(£pl,8B)  in  the  range  of  directions  for  which  -90  deg  <  8  < 

90  deg  or  for  which  index  m  is  in  the  range  46  <  m  <  136  .  Reflected  peak 
direction  is  the  maximum  of  S(£pr,8m)  for  directions  satisfying  8  <  -90  deg 
and  8  >  90  deg  or  direction  indices  in  the  ranges  1  <  m  <  46  and 
136  <  m  <  180.  Note  that  all  reflected  directions  are  searched  for  , 
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and  not  just  those  outside  the  120-deg  azimuths  that  were  used  to  search  for 
fp  r  .  This  avoids  biasing  0p  r  ,  but  does  allow  peaks  occasionally  to  exist 
at  ±90  deg  due  to  the  tails  of  incident  directional  distributions,  as  dis¬ 
cussed  above. 

45.  Reflection  coefficient.  Following  notation  and  definitions  given 
in  the  SPM,  reflection  coefficients  x  ate  generally  defined  as  ratios  of 
reflected  to  incident  wave  heights .  Coefficients  indicating  ratios  of  energy 
would  then  vary  as  the  squares  of  height-based  coefficients  defined  here. 

Three  reflection  coefficients  are  identified.  The  first  is  a  bulk  parameter 
that  characterizes  total  reflected  to  total  incident  energy.  These  two 
energies  are  the  volumes  under  the  reflected  and  incident  spectra,  respective¬ 
ly,  and,  after  removal  of  estimated  noise,  are  the  sources  of  the  definitions 
°f  H-o,r  and  !!,„  t  .  Because  it  is  a  wave -height  ratio  that  is  used  in  the 
definition,  a  bulk  reflection  coefficient  x«o  is  simply 


X-o 


H-o.r 

H-o.i 


(13) 


The  other  two  reflection  coefficients  are  frequency  specific.  One  is  based  on 
the  square  root  of  the  ratio  of  reflected  to  incident  energy  at  the  peak 
frequency  of  the  incident  spectrum  and  suggests  how  much  of  the  energy  at  the 
frequency  spectral  peak  is  reflected.  Denoted  *p.i  ,  it  is  defined  from  the 
incident  and  reflected  frequency  spectra  given  by  Equations  10  and  11, 
respectively,  and  modified  for  noise  removal,  as 


(14) 


where  the  noise  estimate  S„(fp>i)  is  the  average  of  the  smallest  10  values  of 
S(fp>1,0B)  and  the  product  with  180  deg  yields  the  area  under  the  directional 
distribution  attributed  to  noise.  Because  the  peak  in  the  reflected  energy 
may  not  occur  at  the  peak  frequency  of  the  incident  spectrum,  a  third  reflec¬ 
tion  coefficient  is  defined  from  the  ratio  of  frequency  spectral  energies  at 
the  peak  frequency  of  the  reflected  spectrum,  with  a  correction  for  estimated 
noise.  Denoted  xp,r  •  ^  *s  defined  by 
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Sr<f„,r)  -  S,(fp,r)»180- 
Si<fp.r)  -  VW180’ 


(15) 


where  S„(fp  r)  is  the  average  of  the  10  smallest  values  of  S(fp r,0„)  . 

46.  Together,  these  parameters  give  a  bulk  characterization  of  some 
properties  of  a  frequency-direction  spectrum.  There  are,  of  course,  many 
other  parameters  that  can  be  defined,  but  the  present  set  is  simple  and  is 
adequate  for  this  preliminary  analysis. 


Archived  Results 

47.  A  magnetic  storage  medium  containing  the  set  of  observed  frequency - 
direction  spectra  has  been  prepared  for  copying  and  distribution.  These 
spectra  are  immediately  available  to  researchers  within  the  US  Army  Corps  of 
Engineers  and  are  available  to  the  general  public  after  June  1994.  They  can 
be  obtained  by  communicating  with: 

Chief,  Field  Research  Facility 
1261  Duck  Road 
Kitty  Hawk,  NC  27949-4471 
Phone:  (919)  261-3511 

Fax:  (919)  261-4432 

48.  Appendix  A  contains  a  listing  of  the  date,  starting  time  (EST),  and 
the  characterizing  parameters  defined  previously  for  each  case.  It  is 
intended  to  be  used  as  a  kind  of  index  or  catalog  of  the  set  of  cases.  Dates 
are  given  in  the  form  yymmdd  ,  where  yy  is  a  two-digit  year  indicator 
(e.g.,  90  means  1990),  mm  is  the  numeric  index  of  the  calendar  month  (i.e, 

01  is  January,  12  is  December,  etc.),  and  dd  is  the  day  of  the  month.  A 
24-hr  clock  is  used. 

49.  Some  of  the  bulk  parameters  defined  above  are  plotted  as  functions 

of  time  in  Appendix  £.  One  graph  is  shown  for  each  month  of  the  post- 
experiment  period.  Shown  are  H,,,  ,  H,,,^  ,  ,  depth,  tfp  t  ,  a  modified 

form  of  tfp>t  ,  Tp  t  ,  Tp>r  ,  Xmo  .  XPtl  ,  *nd  xp>r  •  The  modification  to  6p  t 

is  a  shift  of  180  deg  so  that,  in  the  figures,  this  variable  represents 

direction  toward  which  reflected  waves  propagate.  The  variable  tfp  t  contin¬ 
ues  to  represent  direction  from  which  incident  waves  arrive. 
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Discuss Ign 


50.  The  intent  of  this  report  is  to  indicate  that  additional  frequency- 
direction  spectral  observations  have  been  added  to  the  already  extensive  FRF 
database  and  to  characterize  these  observations  with  some  simple  parameters. 
There  is  clearly  much  more  work  to  be  done  with  this  extraordinary  data  set, 
including  scientific  interpretation,  but  that  is  beyond  the  scope  of  this 
document.  The  descriptive  parameters  used  here  are  defined  simply  to  be 
preliminary  indicators  of  the  nature  of  the  data.  Close  examination  of  the 
parameters  displayed  in  Appendix  B  indicates  that  the  rudimentary  definitions 
used  here  may  be  somewhat  misleading.  Hence,  a  few  warnings  follow. 

51.  In  plots  of  characteristic  wave  height,  it  is  clear  that  most  of 

the  total  wave  energy,  characterized  by  ,  is  in  the  incident  spectrum, 

which  is  characterized  by  H*,,  A  .  The  so-called  reflected  energy,  indicated 
by  H^,  c  ,  is  small,  but  is  very  seldom  zero.  It  is  possible  that  some  energy 
always  reflects  from  the  FRF  beach,  but  a  more  likely  reason  is  that  the  tails 
of  directional  distributions  for  incident  waves  extend  beyond  the  arbitrarily 
assigned  90-deg  limit  for  incident  energy,  as  mentioned  in  the  paragraphs 
describing  peak  frequency.  Consequently,  i^  t  and  all  the  reflection 
coefficients  are  affected  by  this  characteristic.  Reflection  coefficients  are 
legitimate  estimators  of  seaward -propagating  energy  when  defined  in  the  simple 
way  used  here,  but  a  more  elaborate  scheme  is  required  to  distinguish  direc¬ 
tionally  broad  incident  spectra  from  true  reflected  spectra. 

52.  An  attempt  has  been  made  to  deal  with  the  low,  persistent  levels  of 
energy,  interpreted  here  as  noise,  in  all  the  directional  spectra.  If  this 
interpretation  is  correct,  and  other  sources  of  noise  are  not  dominant,  an 
approximate  confidence  level  can  be  assigned  to  the  reflection  coefficients 
reported  here.  Appendix  C  describes  an  analysis  using  artificial  directional 
distributions  that  include  directionally  broad  noise  at  various  levels,  as 
well  as  Incident  spectra  with  tails  that  cross  the  90-deg  azimuth.  The 
results  of  Appendix  C  suggest  that,  for  noise-to-signal  (N/S)  ratios  estimated 
from  the  present  data  set  to  be  in  the  range  0.01  to  0.1,  moderate  to  high 
estimated  reflection  coefficients  can  be  low  by  as  much  as  25  percent  and  high 
by  as  much  as  50  percent,  with  more  probable  errors  in  the  range  of  5  percent 
to  20  percent  high  for  high-wave  (small  N/S)  conditions.  Small  reflection 
coefficients  (order  0.1)  reported  for  low-wave  conditions  can  be  in  error  by  a 
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factor  of  2  or  more.  Because  the  results  of  Appendix  C  are  based  on  rather 
clean  synthetic  data  and  because  there  nay  be  other  types  of  noise  in  real 
data,  reflection  coefficients  reported  here  must  be  considered  rough  first 
approximations . 

53.  These  observations  qualify,  but  do  not  preclude,  the  usefulness  of 
the  parameters  defined  here.  High  energy  conditions  are  readily  identified, 
as  are  primary  directions  and  periods.  In  high-energy  conditions,  reflected 
wave  heights  rise  sufficiently  above  the  background  level  to  suggest  the 
existence  of  legitimately  quantified  reflected  energy.  For  the  data  set  as  a 
whole,  however,  more  elegant  definitions  clearly  are  required  for  sound 
scientific  interpretation,  a  task  reserved  for  future  work. 

Summary 

54.  Data  from  the  post -experiment  period  of  the  SAMSON  project  have 
provided  an  opportunity  to  add  to  the  FRF  database  of  high- resolution, 
directional -spectral  observations.  Processed  data  have  been  put  in  a  form 
that  will  be  highly  accessible  to  researchers  Interested  in  nearshore  process 
es.  Directional  gage,  directional  analysis  algorithms,  and  definitions  of 
characterizing  parameters  are  described.  The  parameters  used  here  are  first 
approximations,  but  adequate  for  case  distinction  in  this  preliminary  analy¬ 
sis.  More  refined  parametric  definitions  will  be  necessary  for  proper 
scientific  interpretation  of  the  data. 
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18.43 

0.300 

0.331 

0.409 

13.30 

901127 

1900 

0.42 

0.38 

0.11 

2.0  -178.0 

15.63 

18.43 

0.299 

0.283 

0.431 

12.77 

901127 

2200 

0.47 

0.43 

0.12 

-2.0  -176.0 

13.63 

13.63 

0.271 

0.330 

0.330 

13.01 

901128 

0100 

0.43 

0.41 

0.12 

-42.0  -170.0 

13.36 

13.63 

0.290 

0.345 

0.472 

13.65 

901128 

0400 

0.44 

0.41 

0.11 

-38.0  -172.0 

13.36 

15.63 

0.266 

0.267 

0.447 

13.52 

901128 

0700 

0.43 

0.42 

0.11 

6.0  -132.0 

13.56 

13.56 

0.249 

0.343 

0.343 

12.95 

901128 

1000 

0.33 

0.51 

0.12 

8.0  -166.0 

3.24 

13.56 

0.234 

0.133 

0.395 

12.97 

901128 

1300 

0.34 

0.31 

0.12 

2.0  -178.0 

5.24 

15.63 

0.246 

0.231 

0.446 

13.47 

901128 

1600 

0.38 

0.53 

0.12 

-70.0  -164.0 

5.32 

15.63 

0.212 

0.132 

0.376 

13.42 

901128 

1900 

0.48 

0.43 

0.11 

-8.0  -168.0 

13.36 

15.63 

0.234 

0.289 

0.466 

12.76 

901128 

2200 

0.46 

0.43 

0.11 

10.0  -148.0 

13.36 

13.56 

0.250 

0.341 

0.341 

12.76 

901129 

0100 

0.33 

0.31 

0.11 

-34.0  -172.0 

3.24 

13.56 

0.207 

0.173 

0.317 

13.48 

901129 

0400 

0.36 

0.53 

0.13 

-44.0  -172.0 

5.24 

13.36 

0.237 

0.249 

0.238 

13.63 

901129 

0700 

0.32 

0.50 

0.10 

-56.0  -142.0 

5.83 

13.56 

0.205 

0.129 

0.323 

12.97 

901129 

1000 

1.11 

1.07 

0.27 

36.0  90.0 

3.24 

4.75 

0.247 

0.101 

0.289 

12.74 

901129 

1300 

1.69 

1.66 

0.27 

48.0  90.0 

6.19 

4.35 

0.163 

0.086 

0.315 

13.39 

901129 

1600 

1.67 

1.61 

0.37 

36.0  90.0 

5.83 

4.17 

0.230 

0.174 

0.348 

13.70 

901129 

1900 

1.76 

1.70 

0.41 

34.0  90.0 

5.83 

4.17 

0.244 

0.169 

0.393 

13.21 

901129 

2200 

1.93 

1.88 

0.46 

36.0  90.0 

6.19 

5.83 

0.247 

0.221 

0.329 

12.97 

901130 

0100 

2.38 

2.30 

0.35 

26.0  90.0 

7.04 

4.17 

0.241 

0.119 

0.442 

13.71 

901130 

0400 

2.13 

2.03 

0.30 

34.0  90.0 

7.04 

4.17 

0.246 

0.134 

0.383 

14.24 

901130 

0700 

1.86 

1.80 

0.41 

38.0  90.0 

6.39 

4.17 

0.229 

0.164 

0.335 

13.73 

901130 

1000 

1.68 

1.64 

0.27 

30.0  90.0 

7.04 

4.17 

0.167 

0.098 

0.287 

13.03 

901130 

1300 

1.47 

1.44 

0.22 

34.0  90.0 

7.04 

4.17 

0.154 

0.088 

0.272 

13.31 

901130 

1600 

1.30 

1.27 

0.20 

32.0  106.0 

7.04 

7.36 

0.138 

0.134 

0.196 

13.83 

901130 

1900 

1.14 

1.11 

0.16 

16.0  90.0 

8.16 

4.98 

0.143 

0.093 

0.207 

13.40 

901130 

2200 

1.01 

0.99 

0.11 

14.0  90.0 

8.16 

7.56 

0.113 

0.095 

0.131 

12.78 

901201 

0100 

0.88 

0.86 

0.13 

14.0  132.0 

7.04 

8.87 

0.147 

0.108 

0.236 

13.30 

901201 

0400 

0.76 

0.73 

0.13 

22.0  -168.0 

6.19 

7.36 

0.181 

0.131 

0.220 

14.13 

901201 

0700 

0.63 

0.59 

0.13 

14.0  138.0 

7.36 

8.87 

0.213 

0.219 

0.288 

13.79 

901201 

1300 

0.51 

0.49 

0.08 

-28.0  160.0 

9.71 

10.72 

0.167 

0.142 

0.401 

12.82 

901201 

1900 

0.46 

0.43 

0.10 

-28.0  -166.0 

9.71 

11.98 

0.233 

0.143 

0.292 

13.39 

901202 

0117 

0.48 

0.43 

0.11 

10.0  -172.0 

1C.  72 

13.36 

0.244 

C .  22* 

C.108 

12.92 

901202 

0700 

0.41 

0.37 

0.11 

-26.0  -170.0 

9.71 

11.98 

0.293 

0.170 

0.376 

13.92 

901202 

1300 

0.39 

0.37 

0.08 

-32.0  -162.0 

9.71 

11.98 

0.224 

0.134 

0.310 

12.61 

901202 

1900 

0.41 

0.38 

0.10 

-32.0  -160.0 

9.71 

10.72 

0.234 

0.245 

0.243 

13.66 

901203 

0117 

0.46 

0.42 

0.09 

-8.0  -168.0 

9.71 

11.98 

0.222 

0.126 

0.323 

12.89 

901203 

0700 

1.02 

0.98 

0.09 

12.0  170.0 

4.98 

10.72 

0.090 

0.100 

0.193 

14.28 

901203 

1300 

1.10 

1.08 

0.10 

0.0  -90.0 

5.83 

4.98 

0.090 

0.067 

0.116 

12.79 

901203 

1900 

1.37 

1.34 

0.16 

-34.0  -90.0 

7.36 

3.32 

0.119 

0.078 

0.144 

13.86 
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90120* 

0117 

1.56 

1.51 

0.28 

-32.0  -17*. 0 

7.36 

13.63 

0.183 

0.062 

0.13* 

12.67 

90120* 

0700 

1.62 

1.36 

0.32 

-38.0  -178.0 

8.87 

15.63 

0.208 

0.12* 

0.231 

14.07 

90120* 

1300 

1.09 

1.05 

0.16 

-2.0  -90.0 

13.36 

6.39 

0.153 

0.109 

0.305 

12.78 

90120* 

1900 

0.70 

0.7* 

0.1* 

2.0  -150.0 

13.36 

10.72 

0.193 

0.208 

0.263 

13.78 

901205 

0117 

0.S5 

0.82 

0.12 

6.0  -180.0 

11.98 

13.56 

0.150 

0.15* 

0.21* 

12.8* 

901203 

0700 

0.87 

0.83 

0.12 

32.0  -160.0 

6.19 

11.98 

0.133 

0.108 

0.156 

13.99 

901203 

1300 

0.80 

0.79 

0.11 

*0.0  -168.0 

6.19 

11.98 

0.1*0 

0.079 

0.17* 

12.98 

901203 

1900 

0.57 

0.55 

0.08 

36.0  172.0 

3.83 

10.72 

0.132 

0.119 

0.222 

13.50 

901200 

0117 

0.*9 

0.*7 

0.08 

-10.0  -172.0 

11.98 

13.63 

0.177 

0.223 

o.*s* 

12.86 

901200 

0700 

0 .  *8 

0 .  *3 

0.09 

6.0  -176.0 

10.72 

11.98 

0.189 

0.249 

0.297 

13.37 

901200 

1300 

0 .  *6 

0.*3 

0.09 

0.0  -17*. 0 

13.36 

13.36 

0.213 

0.263 

0.265 

13.09 

901200 

1900 

0.*1 

0.38 

0.10 

-* .0  -180.0 

13.36 

13.63 

0.276 

0.233 

0.301 

13.12 

901207 

0117 

0.39 

0.37 

0.09 

-*.0  -172.0 

13.36 

13.36 

0.23* 

0.216 

0.216 

13.09 

901207 

0700 

0.36 

0.33 

0.08 

6.0  -166.0 

13.56 

13.36 

0.2*9 

0.261 

0.261 

13.32 

901207 

1300 

0.37 

0.35 

0.08 

6.0  -166.0 

13.36 

13.56 

0.2*5 

0.230 

0.230 

13.49 

901207 

1900 

0.39 

0.37 

0.09 

-12.0  -166.0 

11.98 

11.98 

0.232 

0.211 

0.211 

13.20 

901200 

0117 

0.76 

0.7* 

0.08 

38.0  -15*. 0 

*.35 

11.98 

0.112 

0.111 

0.135 

13.69 

901200 

0700 

1.33 

1.51 

0.2* 

12.0  90.0 

5.2* 

*.17 

0.138 

0.136 

0.262 

13.3* 

901200 

1300 

2.06 

2.01 

0.31 

6.0  90.0 

6.19 

*.17 

0.156 

0.118 

0.3*7 

13.83 

901200 

1900 

2.03 

1.99 

0.22 

8.0  90.0 

6.39 

*.17 

0.111 

0.071 

0.23* 

13.27 

901209 

0117 

2.07 

2.02 

0.28 

2*  .0  90.0 

6.39 

*.17 

0.136 

0.093 

0.289 

13.98 

901209 

0700 

2.08 

2.03 

0.28 

*.0  90.0 

7.36 

*.17 

0.139 

0.031 

0.332 

13.32 

901209 

1300 

1.7* 

1.70 

0.21 

1*.0  -170.0 

8.16 

8.87 

0.123 

0.078 

0.13* 

13.92 

901209 

1900 

1.27 

1.23 

0.11 

*.0  166.0 

8.16 

3.2* 

0.092 

0.078 

0.145 

13.08 

901210 

0117 

0.92 

0.88 

0.18 

*.0  1*8.0 

8.87 

8.16 

0.210 

0.120 

0.333 

13.82 

901210 

0700 

0.72 

0.69 

0.09 

-2.0  136.0 

8.87 

10.72 

0.129 

0.08* 

0.171 

13.00 

901210 

1300 

0.3* 

0.51 

0.11 

-20.0  170.0 

10.72 

10.72 

0.210 

0.218 

0.218 

13.53 

901210 

1900 

0.31 

0.*8 

0.10 

-8.0  -160.0 

9.71 

11.98 

0.212 

0.106 

0.267 

12.80 

901211 

0117 

1.33 

1.33 

0.16 

3*.0  90.0 

3.32 

*.9S 

0.120 

0.077 

0.165 

13.76 

901211 

0700 

1.26 

1.2* 

0.16 

*0.0  90.0 

6.19 

*.17 

0.126 

0.090 

0.237 

13.13 

901211 

1300 

1.19 

1.17 

0.13 

28.0  130.0 

6.39 

7.36 

0.113 

0.091 

0.137 

13.38 

901211 

1900 

0.91 

0.89 

0.09 

10.0  -168.0 

7.0* 

13.56 

0.102 

0.065 

0.180 

13.01 

901212 

0117 

0.72 

0.70 

0.09 

26.0  -166.0 

3.83 

13.36 

0.126 

0.073 

0.2*3 

13.72 

901212 

0700 

0.39 

0.36 

0.08 

-12.0  -162.0 

11.98 

13.36 

0.1*9 

0.176 

0.320 

13.16 

901212 

1300 

0.*9 

Q.*6 

0.09 

-8.0  -170.0 

13.36 

13.56 

0.199 

0.236 

0.236 

13.38 

901212 

1900 

0.*3 

0.*0 

0.08 

-12.0  -176.0 

13.56 

13.56 

0.212 

0.303 

0.303 

13.01 

901213 

0117 

0 .  *5 

0.*2 

0.10 

-16.0  -162.0 

* .  98 

11.98 

0.228 

0.216 

0.335 

13.46 

901213 

1300 

o.*o 

0.37 

0.08 

-68.0  -168.0 

6.19 

10.72 

0.20* 

0.162 

0.*** 

13.07 

901213 

1900 

0.3* 

0.32 

0.06 

-2*.0  -16*. 0 

6.39 

11.98 

0.197 

0.237 

0.338 

13.06 

90121* 

0117 

1.3* 

1.30 

0.28 

3*  .0  90.0 

*.98 

*.5* 

0.215 

0.223 

0.2*0 

13.48 

90121* 

0700 

2.16 

2.12 

0.33 

36.0  90.0 

7.0* 

*.17 

0.15* 

0.038 

0.283 

13.76 

90121* 

1300 

1.3* 

1.32 

0.16 

8.0  90.0 

7.0* 

*.17 

0.106 

0.075 

0.193 

13.30 

90121* 

1900 

1.2* 

1.22 

0.13 

22.0  90.0 

8.16 

3.2* 

0.107 

0.072 

0.173 

13.** 

901213 

0117 

1.19 

1.17 

0.10 

10.0  90.0 

3.32 

3.2* 

0.090 

0.075 

0.1*6 

13.33 

901213 

0700 

1.01 

0.97 

0.10 

6.0  90.0 

3.32 

*.98 

0.107 

0.128 

0.204 

13.66 

901213 

1300 

1.0* 

1.02 

0.09 

2* .0  -90.0 

6.19 

6.19 

0.08* 

0.081 

0.081 

12.98 

901213 

1900 

1.08 

1.03 

0.13 

-38.0  -90.0 

6.39 

*.98 

0.143 

0.117 

0.213 

13.41 
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901216 

0117 

0.82 

0.80 

0.10 

-44.0 

-90.0 

7.04 

5.32 

0.131 

0.103 

0.212 

13.02 

901216 

0700 

0.62 

0.60 

0.12 

-52.0 

-132.0 

7.04 

8.16 

0.201 

0.188 

0.234 

13.38 

901216 

1300 

0.36 

0.14 

0.06 

-54.0 

-128.0 

7.04 

9.71 

0.187 

0.133 

0.293 

12.75 

901216 

1900 

0.48 

0.46 

0.06 

40.0 

90.0 

4.17 

4.17 

0.133 

0.173 

0.173 

13.47 

901217 

0117 

0.32 

0.51 

0.05 

64.0 

-138.0 

4.17 

7.36 

0.091 

0.112 

0.092 

13.03 

901217 

0700 

0.77 

0.73 

0.08 

44.0 

90.0 

5.83 

5.24 

0.110 

0.109 

0.156 

13.78 

901217 

1300 

0.63 

0.62 

0.06 

20.0 

164.0 

3.83 

11.98 

0.091 

0.036 

0.228 

12.75 

901217 

1900 

0.38 

0.36 

0.06 

20.0 

170.0 

3.83 

5.32 

0.114 

0.093 

0.181 

13.45 

901218 

0117 

0.33 

0.51 

0.04 

20.0 

-170.0 

6.39 

5.83 

0.077 

0.067 

0.130 

12.83 

901218 

0700 

0.46 

0.41 

0.07 

-20.0 

-146.0 

7.36 

8.16 

0.139 

0.172 

0.224 

13.75 

901218 

1300 

0.40 

J .  38 

0.06 

-18.0 

-124.0 

8.16 

11.98 

0.169 

0.149 

0.372 

12.63 

901218 

1900 

0.34 

0.32 

0.12 

-62.0 

-90.0 

4.17 

5.24 

0.232 

0.373 

0.285 

13.30 

901219 

0117 

0.34 

0.32 

0.09 

-48.0 

-90.0 

8.16 

3.52 

0.163 

0.112 

0.270 

12.63 

901219 

0700 

0.33 

0.31 

0.09 

-40.0 

-130.0 

9.71 

9.71 

0.179 

0.161 

0.161 

13.69 

901219 

1300 

0.30 

0.49 

0.05 

-32.0 

178.0 

9.71 

11.98 

0.102 

0.070 

0.319 

12.86 

901219 

1900 

1.03 

1.02 

0.10 

2.0 

-146.0 

4.98 

8.87 

0.100 

0.098 

0.238 

13.64 

901220 

0117 

1.17 

1.13 

0.14 

14.0 

90.0 

5.32 

4.17 

0.121 

0.111 

0.168 

13.02 

901220 

0700 

1.45 

1.43 

0.14 

24.0 

90.0 

5.52 

4.73 

0.098 

0.067 

0.176 

13.95 

901220 

1300 

1.32 

1.30 

0.13 

38.0 

90.0 

6.19 

4.17 

0.087 

0.054 

0.134 

13.13 

901220 

1900 

1.37 

1.55 

0.11 

16.0 

148.0 

6.39 

7.36 

0.06S 

0.034 

0.116 

13.67 

901221 

0117 

1.33 

1.32 

0.10 

16.0 

142.0 

7.04 

8.16 

0.068 

0.042 

0.079 

13.09 

901221 

0700 

1.52 

1.48 

0.12 

-8.0 

-140.0 

7.04 

8.16 

0.082 

0.067 

0.096 

13.79 

901221 

1300 

1.34 

1.50 

0.11 

-12.0 

-90.0 

8.16 

5.32 

0.071 

0.047 

0.136 

13.13 

901221 

1900 

1.36 

1.33 

0.14 

-8.0 

-90.0 

7.36 

5.52 

0.089 

0.067 

0.140 

13.48 

901222 

0117 

1.33 

1.32 

0.13 

-26.0 

-90.0 

8.16 

5.83 

0.099 

0.033 

0.159 

13.14 

901222 

0700 

1.34 

1.31 

0.15 

-8.0 

-130.0 

8.16 

8.16 

0.113 

0.112 

0.112 

13.34 

901222 

1300 

1.05 

1.02 

0.13 

-8.0 

-164.0 

8.16 

11.98 

0.126 

0.078 

0.172 

13.17 

901222 

1900 

1.03 

1.01 

0.12 

-6.0 

-106.0 

7.04 

8.16 

0.116 

0.069 

0.108 

13.30 

901223 

0117 

0.89 

0.86 

0.12 

-8.0 

-146.0 

9.71 

10.72 

0.143 

0.101 

0.143 

13.19 

901223 

0700 

1.10 

1.07 

0.14 

-38.0 

-90.0 

6.39 

5.24 

0.131 

0.075 

0.233 

13.21 

901223 

1300 

1.01 

0.98 

0.14 

-42.0 

-120.0 

7.04 

11.98 

0.149 

0.089 

0.231 

13.14 

901223 

1900 

0.91 

0.88 

0.12 

-14.0 

-116.0 

10.72 

10.72 

0.133 

0.133 

0.133 

12.96 

901224 

0117 

0.81 

0.79 

0.12 

-26.0 

-90.0 

10.72 

3.83 

0.137 

0.113 

0.367 

13.18 

901224 

0700 

0.90 

0.87 

0.14 

-34.0 

-90.0 

10.72 

3.83 

0.161 

0.114 

0.196 

12.91 

901224 

1300 

0.91 

0.87 

0.14 

-16.0 

-132.0 

10.72 

10.72 

0.157 

0.133 

0.153 

13.23 

901224 

1900 

1.84 

1.82 

0.25 

38.0 

90.0 

6.39 

4.17 

0.139 

0.053 

0.293 

13.14 

901223 

.100 

1.35 

1.51 

0.24 

30.0 

90.0 

6.39 

4.17 

0.159 

0.070 

0.330 

13.69 

901223 

0700 

1.17 

1.15 

0.15 

50.0 

90.0 

5.32 

4.17 

0.132 

0.075 

0.236 

13.03 

901223 

1300 

1.09 

1.06 

0.13 

-14.0 

138.0 

10.72 

10.72 

0.122 

0.108 

0.108 

13.43 

901223 

1900 

1.05 

1.02 

0.11 

-16.0 

-90.0 

11.98 

11.98 

0.104 

0.082 

0.082 

12.83 

901226 

0117 

1.02 

0.98 

0.14 

-22.0 

-160.0 

11.98 

13.36 

0.141 

0.118 

0.247 

13.72 

901226 

0700 

0.91 

0.90 

0.11 

-8.0 

-90.0 

11.98 

11.98 

0.117 

0.095 

0.095 

13.03 

901226 

1300 

1.16 

1.12 

0.13 

-6.0 

90.0 

11.98 

4.54 

0.120 

0.110 

0.169 

13.65 

901226 

1900 

1.21 

1.17 

0.13 

16.0 

90.0 

5.24 

4.73 

0.107 

0.094 

0.159 

12.94 

101227 

0117 

1.37 

1.14 

0.14 

12.0 

90.0 

5.24 

4.75 

0.106 

0.095 

0.154 

13.94 

901227 

0700 

1.65 

1.61 

0.24 

16.0 

100.0 

3.52 

4.73 

0.148 

0.137 

0.305 

13.13 

901227 

1300 

1.91 

1.87 

0.23 

8.0 

90.0 

6.19 

6.19 

0.123 

0.162 

0.162 

13.70 

901227 

1900 

1.79 

1.77 

0.12 

14.0 

90.0 

6.39 

4.98 

0.066 

0.046 

0.143 

13.02 
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901228 

0100 

1.88 

1.83 

0.16 

4.0 

-144.0 

7.36 

8.16 

0.088 

0.078 

0.078 

13.94 

901228 

0700 

1.60 

1.36 

0.12 

4.0 

-166.0 

8.87 

8.16 

0.076 

0.046 

0.067 

13.30 

901228 

1300 

1.29 

1.24 

0.15 

-8.0 

-174.0 

8.87 

8.16 

0.118 

0.085 

0.121 

13.64 

901228 

1900 

1.01 

0.97 

0.09 

6.0 

130.0 

8.87 

8.16 

0.094 

0.046 

0.110 

13.13 

901229 

0117 

1.05 

1.00 

0.13 

-2.0 

-172.0 

9.71 

8.16 

0.134 

0.137 

0.132 

13.87 

901229 

0700 

1.07 

1.02 

0.14 

4.0 

160.0 

9.71 

10.72 

0.141 

0.092 

0.209 

13.57 

901229 

1300 

1.13 

1.11 

0.10 

2.0 

140.0 

10.72 

11.98 

0.087 

0.104 

0.133 

13.39 

901229 

1900 

1.05 

1.02 

0.11 

16.0 

-144.0 

10.72 

10.72 

0.109 

0.111 

0.111 

13.28 

901230 

0117 

1.00 

0.97 

0.15 

4.0 

164.0 

10.72 

11.98 

0  134 

0.103 

0.207 

13.43 

901230 

0700 

0.88 

0.83 

0.16 

14.0 

172.0 

10.72 

11.98 

0.196 

0.171 

0.183 

13.61 

901230 

1300 

0.82 

0.79 

0.10 

4.0 

-128.0 

10.72 

9.71 

0.132 

0.079 

0.118 

12.82 

901230 

1900 

0.75 

0.72 

0.15 

-20.0 

-124.0 

10.72 

11.98 

0.203 

0.130 

0.246 

13.19 

901231 

0100 

0.66 

0.62 

0.10 

-8.0 

-138.0 

10.72 

10.72 

0.163 

0.178 

0.178 

12.68 

901231 

1300 

1.65 

1.62 

0.24 

46.0 

90.0 

6.19 

4.17 

0.146 

0.096 

0.249 

12.79 

901231 

1900 

1.44 

1.41 

0.20 

18.0 

90.0 

6.39 

4.17 

0.138 

0.061 

0.233 

13.76 

910101 

0100 

1.38 

1.35 

0.19 

48.0 

90.0 

5.32 

4.33 

0.144 

6.073 

0.223 

12.88 

910101 

0700 

1.40 

1.36 

0.21 

18.0 

-170.0 

3.83 

11.98 

0.134 

0.109 

0.144 

14.23 

910101 

1300 

1.15 

1.13 

0.13 

12.0 

134.0 

6.19 

11.98 

0.111 

0.071 

0.118 

12.84 

910101 

1900 

1.17 

1.13 

0.10 

-8.0 

148.0 

3.32 

11.98 

0.089 

0.091 

0.229 

13.92 

910102 

0117 

1.03 

1.00 

0.10 

-24.0 

144.0 

10.72 

4.33 

0.098 

0.087 

0.144 

12.82 

910102 

0700 

1.15 

1.10 

0.14 

-14.0 

-160.0 

6.19 

11.98 

0.131 

0.131 

0.331 

14.23 

910102 

1300 

1.15 

1.12 

0.15 

24.0 

90.0 

4.35 

4.17 

0.131 

0.226 

0.232 

12.86 

910102 

1900 

1.26 

1.22 

0.11 

26.0 

-160.0 

5.83 

10.72 

0.093 

0.087 

0.267 

13.94 

910103 

0117 

0.92 

0.91 

0.07 

34.0 

-114.0 

3.83 

8.87 

0.079 

0.068 

0.103 

12.80 

910103 

0700 

0.78 

0.75 

0.12 

18.0 

-162.0 

3.83 

8.16 

0.163 

0.140 

0.279 

14.05 

910103 

1300 

0.63 

0.61 

0.08 

-2.0 

138.0 

8.16 

•  16 

0.127 

0.120 

0.120 

12.89 

910103 

1900 

0.72 

0.69 

0.11 

•38.0 

-140.0 

9.71 

8.87 

0.165 

0.198 

0.226 

13.78 

910104 

0100 

0.37 

0.34 

0.09 

-34.0 

162.0 

9.71 

10.72 

0.162 

0.087 

0.227 

13.04 

910104 

0700 

1.20 

1.18 

0.13 

28.0 

90.0 

4.98 

4.17 

0.129 

0.138 

0.226 

13.95 

910104 

1300 

1.34 

1.31 

0.17 

26.0 

90.0 

3.83 

4.35 

0.130 

0.068 

0.199 

13.23 

910104 

1900 

1.05 

1.03 

0.10 

8.0 

90.0 

3.32 

4.35 

0.093 

0.063 

0.133 

13.64 

910103 

0117 

1.06 

1.04 

0.10 

12.0 

90.0 

3.32 

4.17 

0.097 

0.088 

0.102 

13.30 

910103 

0700 

1.17 

1.13 

0.13 

4.0 

90.0 

4.98 

3.24 

0.109 

0.141 

0.127 

13.62 

910103 

1300 

1.02 

1.00 

0.09 

0.0 

90.0 

3.32 

4.98 

0.093 

0.063 

0.144 

13.37 

910103 

1900 

0.96 

0.93 

0.08 

-14.0 

-170.0 

6.39 

8.87 

0.089 

0.072 

0.133 

13.36 

910106 

0117 

0.73 

0.70 

0.11 

-22.0 

-124.0 

6.39 

8.16 

0.160 

0.131 

0.256 

13.46 

910106 

0700 

0.61 

0.39 

0.08 

-4.0 

-90.0 

6.19 

3.52 

0.132 

0.117 

0.133 

13.22 

910106 

1300 

0.47 

0.44 

0.10 

-34.0 

-90.0 

9.71 

5.32 

0.218 

0.149 

0.347 

13.40 

910106 

1900 

0.30 

0.47 

0.07 

-36.0 

-178.0 

8.87 

18.43 

0.144 

0.136 

0.382 

13.06 

910107 

0117 

0.31 

0.48 

0.10 

-38.0 

-90.0 

9.71 

8.87 

0.217 

0.181 

0.272 

13.63 

910107 

0700 

0.30 

0.47 

0.09 

-40.0 

-180.0 

9.71 

18.43 

0.188 

0.172 

0.707 

13.22 

910107 

1300 

1.96 

1.91 

0.30 

34.0 

90.0 

6.59 

4.54 

0.137 

0.089 

0.241 

13.83 

910107 

1900 

2.89 

2.85 

0.34 

8.0 

90.0 

8.16 

4.17 

0.120 

0.083 

0.262 

13.44 

910108 

0117 

2.99 

2.93 

0.41 

10.0 

90.0 

8.16 

4.17 

0.139 

0.084 

0.339 

14.17 

910108 

0700 

2.73 

2.68 

0.29 

0.0 

90.0 

8.87 

4.17 

0.109 

0.066 

0.272 

13.47 

910108 

1300 

3.00 

2.96 

0.23 

4.0 

130.0 

8.16 

4.17 

0.083 

0.034 

0.238 

13.98 

910108 

1900 

2.81 

2.76 

0.25 

4.0 

132.0 

10.72 

4.17 

0.090 

0.071 

0.311 

13.36 
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910109 

0117 

3.03 

2.97 

0.26 

2.0  132.0 

10.72 

4.17 

0.080 

0.044 

0.399 

14.17 

910109 

0700 

2.09 

2.03 

0.31 

-6.0  90.0 

11.98 

4.17 

0.118 

0.046 

0.470 

13.43 

910109 

1300 

2.59 

2.53 

0.26 

-2.0  146.0 

11.98 

10.72 

0.105 

0.042 

0.092 

13.88 

910109 

1900 

1.03 

1.70 

0.14 

-2.0  142.0 

11.98 

11.98 

0.077 

0.070 

0.070 

13.32 

910110 

0117 

1.06 

1.00 

0.21 

6.0  160.0 

10.72 

10.72 

0.132 

0.128 

0.128 

14.03 

910110 

0700 

1.31 

1.46 

0.16 

-20.0  174.0 

9.71 

10.72 

0.112 

0.071 

0.134 

13.40 

910110 

1300 

1.49 

1.45 

0.10 

10.0  -162.0 

10.72 

10.72 

0.123 

0.175 

0.173 

13.73 

910110 

1900 

1.13 

1.11 

0.13 

8.0  -134.0 

10.72 

10.72 

0.116 

0.172 

0.172 

13.29 

910111 

0117 

1.49 

1.43 

0.17 

4.0  -172.0 

3.24 

10.72 

0.120 

0.119 

0.164 

13.92 

910111 

0700 

1.62 

1.50 

0.15 

20.0  -178.0 

0.39 

10.72 

0.093 

0.067 

0.192 

13.52 

910111 

1300 

1.30 

1.52 

0.16 

8.0  -170.0 

7.04 

9.71 

0.100 

0.000 

0.105 

13.50 

910111 

1900 

1.99 

1.93 

0.20 

-6.0  -90.0 

8.10 

8.87 

0.105 

0.079 

0.077 

13.28 

910112 

0117 

2.33 

2.30 

0.20 

-32.0  -90.0 

8.87 

8.16 

0.111 

0.000 

0.110 

13.60 

910112 

0700 

1.00 

1.00 

0.22 

-30.0  176.0 

9.71 

10.72 

0.140 

0.097 

0.174 

13.33 

910112 

1300 

1.40 

1.35 

0.20 

-32.0  -116.0 

8.87 

8.16 

0.149 

0.088 

0.150 

13.33 

910112 

1900 

1.13 

1.09 

0.21 

-32.0  -164.0 

8.87 

10.72 

0.189 

0.112 

0.200 

13.41 

910113 

0117 

1.21 

1.13 

0.14 

-32.0  90.0 

8.87 

4.17 

0.120 

0.102 

0.269 

13.38 

910113 

0700 

1.42 

1.36 

0.26 

28.0  90.0 

3.32 

4.33 

0.192 

0.154 

0.299 

13.81 

910113 

1300 

1.13 

1.09 

0.13 

32.0  90.0 

5.83 

4.75 

0.117 

0.073 

0.201 

13.29 

910113 

1900 

1.07 

1.03 

0.10 

34.0  132.0 

5.83 

4.54 

0.179 

0.132 

0.294 

13.38 

91011* 

0117 

0.07 

0.04 

0.11 

30.0  90.0 

5.83 

4.75 

0.133 

0.077 

0.229 

13.44 

910114 

0700 

0.74 

0.09 

0.12 

20.0  176.0 

5.83 

3.24 

0.181 

0.134 

0.206 

13.87 

910114 

1300 

0.33 

0.52 

0.00 

-28.0  -178.0 

9.71 

15.63 

0.132 

0.090 

0.342 

13.04 

910114 

1900 

0.33 

0.30 

0.11 

-0.0  -134.0 

13.03 

18.43 

0.223 

0.265 

0.306 

13.40 

910113 

0117 

0.49 

0.46 

0.09 

-6.0  -178.0 

15.03 

15.03 

0.200 

0.230 

0.236 

13.01 

910113 

0700 

0.45 

0.41 

0.11 

-8.0  -174.0 

13.03 

15.63 

0.233 

0.236 

0.230 

13.73 

910113 

1300 

0.46 

0.41 

0.11 

-8.0  -172.0 

15.03 

15.03 

0.230 

0.375 

0.375 

12.80 

910113 

1900 

0.43 

0.42 

0.10 

6.0  -180.0 

13.50 

15.03 

0.243 

0.239 

0.388 

13.00 

910110 

0117 

0.43 

0.42 

0.00 

-6.0  -174.0 

13.36 

13.56 

0.202 

0.233 

0.235 

13.02 

910116 

0700 

1.74 

1.67 

0.33 

-20.0  -126.0 

8.10 

8.16 

0.200 

0.152 

0.152 

13.80 

910110 

1300 

1.42 

1.30 

0.17 

-34.0  -90.0 

8.87 

8.16 

0.124 

0.079 

0.133 

12.83 

910116 

1900 

1.07 

1.02 

0.23 

-50.0  -90.0 

8.10 

8.10 

0.223 

0.210 

0.216 

13.08 

910117 

0117 

0.70 

0.72 

0.14 

-38.0  -140.0 

8.87 

10.72 

0.200 

0.142 

0.321 

12.92 

910117 

0700 

0.63 

0.01 

0.17 

-40.0  -148.0 

9.71 

10.72 

0.284 

0.218 

0.353 

13.90 

910117 

1300 

0.43 

0.42 

0.00 

•38.0  -132.0 

8.87 

10.72 

0.187 

0.111 

0.301 

12.84 

910117 

1900 

0.43 

0.40 

0.09 

-32.0  -90.0 

8.87 

8.87 

0.237 

0.236 

0.230 

13.70 

910110 

0117 

0.36 

0.33 

0.06 

-34.0  168.0 

8.87 

10.72 

0.194 

0.116 

0.323 

12.83 

910110 

0700 

0.37 

0.34 

0.06 

-34.0  -122.0 

9.71 

8.16 

0.180 

0.171 

0.288 

13.73 

910110 

1300 

0.30 

0.27 

0.03 

8.0  170.0 

10.72 

10.72 

0.184 

0.279 

0.279 

12.71 

910110 

1900 

0.23 

0.21 

0.06 

-2.0  -178.0 

9.71 

18.45 

0.282 

0.414 

0.715 

13.34 

910119 

0117 

0.10 

0.10 

0.04 

-32.0  -176.0 

9.71 

18.43 

0.249 

0.233 

0.512 

12.70 

910119 

0700 

0.21 

0.10 

0.06 

28.0  -176.0 

0.19 

18.43 

0.302 

0.135 

0.733 

13.00 

910119 

1300 

0.23 

0.23 

0.04 

-34.0  -180.0 

9.71 

13.63 

0.193 

0.183 

0.518 

12.73 

910119 

1900 

0.20 

0.23 

0.06 

14.0  -170.0 

0.19 

18.45 

0.239 

0.123 

0.072 

13.37 

910120 

0117 

0.23 

0.23 

0.05 

0.0  -176.0 

13.03 

18.45 

0.236 

0.374 

0.013 

12.80 

910120 

0700 

0.47 

0.43 

0.11 

32.0  -180.0 

4.17 

13.03 

0.231 

0.222 

0.055 

13.53 

910120 

1300 

0.70 

0.73 

0.10 

22.0  90.0 

4.34 

4.34 

0.130 

0.100 

0.100 

12.97 

910120 

1900 

0.02 

0.30 

0.14 

-38.0  -128.0 

9.71 

11.98 

0.244 

0.203 

0.352 

13.40 
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910121 

0117 

0.33 

0.30 

0.16 

-20.0 

-130.0 

10.72 

10.72 

0.321 

0.470 

0.470 

13.17 

910121 

0700 

0.46 

0.43 

0.11 

-24.0 

-136.0 

10.72 

9.71 

0.261 

0.242 

0.297 

13.39 

910121 

1300 

0.32 

0.49 

0.12 

60.0 

90.0 

4.35 

4.17 

0.232 

0.339 

0.382 

13.26 

910121 

1900 

l.*4 

1.39 

0.32 

32.0 

90.0 

3.83 

4.98 

0.233 

0.181 

0.307 

13.65 

910122 

0117 

1.33 

1.27 

0.30 

30.0 

90.0 

6.19 

4.17 

0.238 

0.157 

0.319 

13.75 

910122 

0700 

1.39 

1.33 

0.33 

30.0 

106.0 

5.32 

4.98 

0.218 

0.130 

0.356 

13.63 

910122 

1300 

1.6* 

1.39 

0.32 

40.0 

90.0 

6.39 

4.17 

0.204 

0.103 

0.402 

13.69 

910122 

1900 

1.29 

1.23 

0.21 

14.0 

132.0 

7.36 

8.16 

0.168 

0.138 

0.238 

13.4* 

910123 

0117 

0.93 

0.90 

0.15 

32.0 

144.0 

7.04 

7.36 

0.165 

0.130 

0.211 

13.82 
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Appendix  C:  Estimation  of  Error  in  Reflection  Coefficients 


1.  To  obtain  some  measure  of  reliability  for  the  reflection  coeffi¬ 
cients  deduced  from  the  Sources  of  Ambient  Micro-Seismic  Noise  (SAMSON)  data, 
tests  have  been  conducted  using  a  variety  of  synthetic  directional  spectra  in 
conjunction  with  the  Iterative  Maximum  Likelihood  Estimation  (IMLE)  algorithm 
and  the  geometry  of  the  SAMSON  array  of  gages.  Artificial  directional  energy 
distributions  are  used  to  compute  the  cross -spectral  matrices  that  would 
result  from  the  SAMSON  gages  if  the  corresponding  idealized  wave  fields  had 
actually  been  measured.  These  cross -spectral  matrices  are  then  passed  through 
the  same  directional  estimation  algorithm  that  was  used  to  process  the  SAMSON 
field  data.  The  resulting  IMLE  estimates  of  directional  distributions  are 
then  compared  to  the  original  synthetic  distributions  to  assess  errors  in  the 
analysis  algorithm. 

2.  To  keep  the  test  results  tractable,  synthetic  distributions  are 
necessarily  kept  simple.  Consequently,  the  test  distributions  do  not  mimic 
exactly  the  structurally  more  complicated  distributions  detected  in  the  field 
data  reported  in  the  main  body  of  this  report.  As  a  result,  these  tests  are 
only  indicative  of  the  viability  of  the  analysis  algorithm  and  are  not 
specifically  applicable  to  any  particular  real  observation.  Nonetheless,  it 
is  of  immense  value  to  have  some  quantitative  measure  of  efficacy  of  the 
analysis  system,  and  describing  such  a  measure  is  the  intent  of  this  appendix. 

3.  Test  spectra  are  single  directional  distributions  constructed  of  two 
Gaussian- shaped  modes,  one  representing  incident  waves  and  the  other  repre¬ 
senting  reflected  waves,  to  which  is  added  a  white  (constant  with  respect  to 
direction)  noise  spectrum.  A  single  directional  distribution  is  equivalent  to 
a  directionwise  slice  at  a  particular  frequency  through  a  full  frequency- 
direction  spectrum,  of  which  Figure  6  in  the  main  text  is  an  example. 

Frequency  is  specified  for  the  tests  because  it  determines  the  wavelengths  of 
the  (synthetic)  waves  that  would  pass  through  the  array  of  gages,  and  it 
determines  which  sub-array  of  gages  is  used  in  analysis,  following  Table  1  in 
the  main  text.  Computations  are  done  at  discrete  directions*  6a  as  in  the 
main  text,  but  with  the  direction  increment  66  -  1  deg  instead  of  2  deg,  and 
increasing  from  -180  deg.  Hence,  directions  are  defined  by 


A  notation  for  this  appendix  is  included  in  Appendix  D. 


em  *  -180*  +  (m  -  1)  d8  , 


1-1,  2,  3,  ...  .  M 


(Cl) 


where  M  -  360  is  the  total  number  of  direction  increments.  In  terms  of 
these  discrete  directions,  synthetic  spectra  take  the  form 


where  the  first  term  on  the  right-hand  side  represents  the  incident  mode,  the 
second  term  on  the  right-hand  side  is  the  reflected  mode,  and  the  last  term  is 
the  (constant)  noise  spectrum.  The  sum  of  the  first  two  terms  on  the  right- 
hand  side  is  the  signal.  There  are  six  variable  parameters  in  Equation  C2 : 

PA  is  the  fraction  of  total  signal  energy  in  the  incident  mode,  8p  t  is  the 
incident  mode  peak  direction,  68i  governs  the  directional  spread  of  the 
incident  mode,  8p  r  is  the  reflected  mode  peak  direction,  68 r  determines  the 
directional  spread  of  the  reflected  mode ,  and  SN  is  the  total  variance  in 
the  noise  spectrum.  The  coefficient  (1  -  Pi)  in  the  second  term  gives  the 
fractional  part  of  total  signal  energy  in  the  reflected  mode.  Total  noise  SH 
is  divided  by  360  in  Equation  C2  to  give  the  proper  variance  density  for  each 
discrete  direction  8a  .  The  parameters  68i  and  68t  are  like  the  standard 
deviations  in  Gaussian  probability  density  functions,  and  determine  the  rate 
at  which  the  Gaussian  curve  decays  from  its  peak.  In  practice,  these  two 
parameters  are  defined  from  the  expressions 
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and 
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and  t8i  and  A 8t  are  specified  because  they  give  the  actual  directional 
spreads  at  1/e  times  the  peak  values  of  the  Gaussian  curves. 

4.  The  total  variance  represented  by  Equation  C2  is  found  by  integra¬ 
tion  of  both  sides  with  respect  to  8  ,  which,  in  discrete  form,  is  equivalent 
to  summing  the  product  S(0„)  d 8  over  all  directions.  The  total  signal 
variance  has  been  arbitrarily  set  to  1  m2.  The  coefficients  have  been  defined 
so  that  this  discrete  summation  o»er  a  finite  range  of  directions  very  nearly 
approximates  the  integral  of  the  continuous  (in  8  )  form  of  Equation  C2  over 
an  infinite  range  of  directions,  i.e.. 
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In  this  expression,  the  ratio  of  the  total  variance  due  to  noise  (  SK  )  to 
total  signal  variance  (1  m2)  is  simply  SH  ,  so  that  the  noise- to- signal  ratio 
N/S  is  SN  identically.  In  computations  with  the  spectrum  of  Equation  C2, 
the  area  of  the  signal  part  of  the  spectrum  is  computed  separately,  and  the 
first  two  terms  on  the  right-hand  side  of  Equation  C2  are  renormalized  by  this 
area  to  make  Equation  C4  exact  (to  computer  accuracy). 

5.  The  formulation  of  Equation  C2  also  allows  some  (approximate) 
control  of  reflection  coefficients  for  the  true  spectra.  If  peak  directions 
and  directional  spreads  are  such  that  the  tails  of  neither  incident  nor 
reflected  mode  are  significant  on  opposing  sides  of  the  ±90-deg  boundary  that 
separates  incident  and  reflected  wave  energy,  the  wave -height -based  reflection 
coefficient  used  here  will  be  simply 
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However,  because  many  of  the  test  spectra  do  have  extensive  tail  regions, 
reflection  coefficients  for  true  and  estimated  spectra  are  computed  numerical¬ 
ly  in  the  same  way  as  described  in  the  main  text,  modified  for  the  smaller 
direction  increment  dd  .  With  that  modification,  and  following  the  notation 
of  Equation  C2,  the  total  incident  energy  plus  noise  is 
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and  the  total  reflected  energy  plus  noise  is 
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The  reflection  coefficient  is  then  computed  as 
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where  half  the  noise  variance  is  associated  with  the  total  incident  energy  and 
the  other  half  is  aligned  with  total  reflected  energy.  For  the  true  spectrum, 
S„  is  specified  exactly,  so  the  true  reflection  coefficient  can  be  computed 
with  high  accuracy.  For  spectra  estimated  with  the  Maximum  Likelihood 
Estimation  (MLE)  or  IMLE  algorithms,  SN  is  estimated  as  the  average  of  the 
smallest  10  discrete  spectral  values,  just  as  was  done  with  the  SAMSON  field 
data.  The  accuracy  of  reflection  coefficients  computed  from  either  of  the 
estimation  algorithms  thus  depends  strongly  on  how  well  they  can  resolve  the 
noise,  especially  at  high  noise  levels. 

6.  Variation  of  the  six  parameters  Pt  ,  i  ,  9f  r  ,  6^  ,  6BZ  ,  and 
SH  allows  considerable  flexibility  in  test  conditions  for  the  estimation 
algorithms.  Once  these  parameters  are  specified,  Equation  C2  can  be  evaluat¬ 
ed.  When  a  frequency  and  water  depth  are  specified  for  the  test  condition, 
the  set  of  gages  in  the  sub-array  can  be  identified  as  in  Table  1  in  the  main 
text,  and  a  wavelength  A  can  be  found  using  the  dispersion  relation  of 
linear  wave  theory.  The  true  complex  cross -spectral  matrix  element 

Cpq  '  iQpq  between  one  gage  located  at  coordinates  (Xp,yp)  and  another  gage 
located  at  coordinates  (xq,yq)  is  computed  from  the  expression 

Cpq  ~  iQpq  *  Y.  S(0J  g-iktCXp-  xq)co*«m-  (yp-y.jj.infl,,]  ^  (C9) 

where  k  —  2x/A  is  the  wave  number  and  i  is  the  square  root  of  -1. 

Equation  C 9  is  a  discrete  approximation  of  the  exact  sea-surface  variance 
cross-spectral  density  derived  from  first  principles  by  Munk  et  al.  (1963)*. 

The  inversion  of  Equation  C9  to  isolate  S(0B)  is  fundamentally  central  to 
many  directional  estimation  algorithms. 

7.  The  complete  cross-spectral  matrix  is  formed  when  the  computation  of 

Equation  C9  is  performed  for  all  pairs  of  gages  in  the  sub-array.  This  matrix 
is  then  treated  as  if  it  was  derived  from  field  data,  and  processed  with  the 
IMLE  algorithm  described  in  the  main  text.  The  result  is  a  discrete  estimate 
■“*1  ( ^o)  °f  the  true  spectrum  of  Equation  C2,  and  various  comparison  tests  can 
then  be  performed.  In  the  computation  of  an  IMLE  estimate,  the  first  itera¬ 
tion  is  an  MLE  estimate,  denoted  .  It  is  useful  to  have  this  estimate 

as  well,  so  it  is  included  in  some  of  the  comparison  tests  described  below. 


See  References  at  the  end  of  the  main  text. 
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8.  An  example  of  the  type  of  results  obtained  with  these  computations 
is  shown  in  Figure  Cl,  where  the  MLE-  and  IMLE- estimated  directional  distribu¬ 
tions  are  plotted  with  the  true  spectrum.  Various  parameters  describe  the 
goodness  of  fit  of  the  two  estimated  curves  and  the  behavior  of  the  IMLE 
analysis  algorithm.  Parameters  of  the  true  spectrum  are  shown  in  the  page 
header.  Mode  1  information  relates  to  the  incident  signal  spectrum  and  mode  2 
information  is  for  the  reflected  signal  spectrum.  The  entry  in  the  column 
labeled  "size"  for  mode  1  is  the  parameter  PA  ,  and  that  for  mode  2  is 

1  -  Pi  .  The  column  labeled  "peak"  lists  0  A  and  8p  r  for  modes  1  and  2, 
respectively.  The  column  labeled  "spread"  lists  A0a  and  A0r  for  modes  1 
and  2,  respectively.  The  value  shown  for  the  N/S  ratio  is  SK  .  Frequency, 
depth,  and  the  corresponding  wavelength  for  this  test  case  are  shown  in  the 
first  line  of  the  page  header. 

9.  The  line  labeled  "gpats"  contains  two  strings  of  letters,  which  are 
codes  for  subsets  of  gages  from  the  SAMSON  array.  Each  letter  represents  the 
gage  number  corresponding  to  that  letter's  ordinal  position  in  the  alphabet. 
Hence,  the  letter  "A"  represents  gage  1,  "B"  represents  gage  2,  and  so  on 
through  letter  "X",  which  is  gage  24.  Gages  12  (letter  ".")  and  14  (letter 
"N")  are  not  used  in  these  tests  because  they  were  not  used  in  processing  the 
SAMSON  field  data.  The  first  string  of  letters  represents  the  initial  sub¬ 
array  of  gages  appropriate  for  the  wave  frequency  of  this  test,  following 
Table  1  in  the  main  text.  The  second  string  of  letters  represents  the  set  of 
gages  actually  used  to  compute  the  estimated  spectra  shown.  This  second  gage 
set  was  identified  by  the  dynamic  gage  elimination  procedure  described  in  the 
main  text.  In  the  course  of  the  test  computations  described  here,  it  was 
commonly  found  that  the  initial  gage  set  was  reduced  substantially  through 
dynamic  elimination.  This  behavior  occurred  because  cross-spectral  informa¬ 
tion  tended  to  peak  with  relatively  few  gages  for  the  simple  shapes  of  test 
spectra  used  in  these  tests.  The  page  header  in  Figure  Cl  indicates  that 
after  the  final  set  of  gages  was  identified,  it  took  25  iterations  of  the  IMLE 
algorithm  to  converge  to  the  result  shown. 

10.  Computed  parameters  and  parametric  comparisons  are  shown  in  the 
header  inside  the  graph  frame.  The  column  marked  "x"  contains  the  true 
reflection  coefficient  x  .  and  the  reflection  coefficients  from  the  MLE 
estimate  Xm  and  the  IMLE  estimate  Xi  •  Errors  in  estimated  reflection 
coefficients  are  denoted  Ax  ,  and  are  computed  as  percentage  differences  from 
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the  true  reflection  coefficient.  Errors  in  the  MLE  estimate  Ax«  and  IMLE 
estimate  Axi  are  computed  from 

Axm  *  — - -  *  100  percent  and  AXj  *  — - -  *  100  percent  (CIO) 

X  X 

where  X  .  Xm  <  and  Xi  are  found  using  Equations  C6,  C7,  and  C8.  While 
reflection  coefficient  estimation  is  the  subject  of  this  appendix,  some  other 
measures  of  estimator  quality  are  also  given  in  Figure  Cl . 

11.  A  measure  of  overall  fidelity  of  an  estimated  spectrum  to  the  true 
spectrum  is  indicated  by  the  square  root  of  the  mean  squared  difference  of  two 
spectra,  normalized  by  the  peak  of  the  true  spectrum,  and  expressed  as  a 
percent.  The  column  marked  "RMS"  shows  the  results  of  such  a  computation.  If 
RMSh  and  RMSX  denote  these  values  for  the  MLE  and  IMLE  estimates,  respec¬ 
tively,  the  computation  for  RMSj  obeys 

**»!  *  ^  ^  £  [  V*.)  -  S(«ri)f  (cll) 

Ssmx  M  »■! 

where  S„„  [  -  S (0Pii)  ,  in  the  tests  reported  here]  is  the  peak  value  of  the 

true  spectrum.  The  computation  for  RMSm  also  obeys  Equation  Cll,  but  with 
SM  replacing  S:  .  Note  that  values  for  RMSm  and  RMSj  can  appear  quite 
small  if  there  are  large  regions  of  the  spectrum  with  very  low  signal,  as 
illustrated  by  Figure  Cl.  Squared  differences,  sometimes  large,  in  the  signal 
parts  of  the  spectra  can  be  averaged  with  large  numbers  of  small  values  where 
there  is  no  signal,  and  the  resulting  computation  can  give  a  small  number. 

Some  caution  is  therefore  required  in  interpreting  this  measure  of  error. 

12.  Because  bimodal  signal  spectra  are  used  in  the  tests  reported  here, 
the  two  estimated  spectra  can  be  searched  for  incident  and  reflected  peak 
directions.  The  ability  to  resolve  peak  directions  can  be  quantified  by 
comparing  estimated  peaks  to  true  peaks.  The  column  marked  "0p  t"  in  Fig¬ 
ure  Cl  shows  the  true  incident  peak  direction  and  the  peak  incident  directions 
of  the  MLE  and  IMLE  estimates.  The  column  marked  "0p  r"  shows  the  correspond¬ 
ing  reflected  peak  directions.  The  true  incident  and  reflected  peak  direc¬ 
tions  are  the  modal  peak  directions  defined  in  the  page  header,  and  can  always 
be  specified.  In  some  test  cases,  the  modes  can  be  combined  in  such  a  way 
that  there  is  no  specific  peak  for  the  reflected  spectrum,  and,  if  one  or  both 
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of  the  estimated  spectra  faithfully  represent  the  true  spectrum,  no  peak 
reflected  direction  can  be  found.  In  such  cases,  an  entry  of  "999"  is  used  in 
the  column  for  peak  reflected  directions. 

13.  Examination  of  some  of  the  test  parameters  shown  in  Figure  Cl 
suggests  that  the  MLE  estimate  may  be  slightly  better  than  the  IMLE  estimate 
for  this  case.  The  RMS  error  is  slightly  smaller  and  differences  between  true 
and  estimated  incident  and  reflected  peak  directions  are  smaller  for  the  MLE 
estimate.  However,  the  reflection  coefficient  estimate  is  much  worse  for  the 
MLE  estimate  than  for  the  IMLE  estimate.  This  is  because  the  IMLE  algorithm 
is  better  at  distributing  the  estimated  energy  in  direction,  i.e.,  it  is 
better  at  getting  the  right  amounts  of  energy  in  the  semicircles  defined  as 
incident  and  reflection.  As  can  be  seen  in  Figure  Cl,  the  MLE  curve  underes¬ 
timates  the  area  in  the  incident  peak  and  overestimates  the  area  in  the 
reflected  peak  more  than  the  IMLE  estimate  does.  Consequently,  the  ratio  of 
these  two  areas  (which  leads  to  the  reflection  coefficient  in  low-noise  cases) 
is  much  better  represented  by  the  IMLE  estimate,  which,  in  the  case  of  Figure 
Cl,  yields  a  reflection  coefficient  only  about  2  percent  different  than  the 
true  one. 

14.  If  the  noise  level  is  raised  for  the  case  depicted  in  Figure  Cl, 
some  effects  of  a  finite,  directionally  white,  background  signal  on  reflection 
coefficient  estimation  can  be  seen.  Figure  C2  shows  a  case  where  the  true 
signal  properties  are  the  same  as  in  Figure  Cl  but  N/S  is  raised  from  0.0001 
to  0.1,  a  level  high  enough  to  be  seen  as  the  nearly  constant  intermodal  line 
at  the  bottom  of  the  figure.  The  IMLE  estimated  spectrum  rather  faithfully 
reproduces  the  whole  true  spectrum,  and  overpredicts  the  reflection  coeffi¬ 
cient  by  only  about  5  percent.  Fart  of  the  reason  for  this  error  is  that  the 
IMLE  estimate  undulates  slightly  about  the  background  noise  level.  The  method 
of  using  the  average  of  the  smallest  10  points  of  the  discrete  estimated 
spectrum  to  estimate  the  noise  is  probably  biased  low  by  the  low  points  in  the 
undulations.  Consequently,  the  true  noise  level  is  underestimated,  the 
difference  makes  a  greater  contribution  to  the  estimate  of  reflected  area,  and 
the  reflection  coefficient  is  slightly  overestimated.  That  the  error  is  small 
is  somewhat  remarkable.  The  area  under  the  noise  spectrum  is  10  percent  of 
the  signal,  and  the  area  under  the  reflected  peak  (0.1)  is  10  percent  of  the 
total  signal  (1.0),  so  the  ratio  of  noise  to  reflected  signal  is  1.0,  which  is 
rather  high.  On  the  other  hand,  the  reflection  coefficient  is  moderately 
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Figure  C2.  Sane  conditions  as  in  Figure  Cl,  but  with  increased  noise 
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high,  given  the  range  of  reflection  coefficients  listed  in  Appendix  A  and 
depicted  in  Appendix  B. 

15.  Many  of  the  values  listed  in  Appendix  A  are  nearer  to  0.1  than  0.3, 
so  the  test  conditions  shown  in  Figure  C2  can  be  modified  by  changing  the 
relative  sizes  of  incident  and  reflected  energies,  keeping  all  other  governing 
parameters  the  same.  Figure  C3  shows  such  a  case  where  the  incident  energy  is 
raised  to  98  percent  of  the  total  signal,  the  reflected  energy  is  2  percent  of 
the  total  signal,  and  the  true  reflection  coefficient  becomes  approximately 
0.14.  The  behavior  of  the  IMLE  estimate  in  Figure  C3  is  largely  the  same  as 
in  Figure  C2,  but  now  the  slight  errors  mentioned  in  the  last  paragraph  have  a 
larger  effect  on  the  estimate  of  the  reflected  energy,  and  the  reflection 
coefficient  is  overpredicted  by  about  22  percent. 

16.  Uhen  the  modal  peaks  are  close  together,  it  is  sometimes  more 

difficult  for  the  IMLE  algorithm  to  resolve  the  true  spectral  structure. 

Figure  C4  shows  a  rather  extreme  case  where  the  modal  spreads  stay  the  same  as 

in  Figure  C3,  but  the  peak  directions  are  moved  to  within  60  deg  of  each 
other.  Also,  the  area  of  the  reflected  mode  is  halved  from  the  value  used  in 
Figure  C3,  reducing  the  true  reflection  coefficient  to  about  0.11.  The  small 
variations  in  the  IMLE  estimate  from  the  true  spectrum  result  In  an  over¬ 
prediction  of  the  reflection  coefficient  by  about  42  percent.  This  result  is 

not  surprising  because,  in  this  case,  the  area  of  the  noise  spectrum  is  10 
times  the  area  of  the  reflected  signal.  This  suggests  that  field  data  having 
estimated  reflection  coefficients  of  order  0.1,  under  low  energy  (high  N/S) 
conditions,  and  having  high  incident  angles  of  attack  will  be  subject  to  this 
much  error,  at  least.  Additional  errors  (e.g.,  those  due  to  uncertainty  in 
cross -spectral  estimates)  in  field  data  and  not  evaluated  here  will  likely 
increase  errors  in  reflection  coefficients. 

17.  The  errors  shown  in  Figure  C4  are  extreme  for  the  set  of  test 
conditions  evaluated  at  this  frequency  and  noise  level.  If  the  modal  peaks 
remain  close  together,  but  the  modal  spreads  increase  so  that  modal  tail 
regions  both  overlap  and  cross  the  nominal  90-deg  boundary  that  distinguishes 
"incident"  from  "reflected"  energy,  the  errors  decrease  again.  Figure  C5 
shows  a  low-noise  (N/S  -  0.0001)  case  with  modal  peaks  separated  by  60  deg, 
modal  spreads  increased  to  40  deg  (from  20  deg  in  Figure  C4)  and  the  reflected 
mode  given  10  percent  of  total  signal  variance.  Both  modal  tail  regions  cross 
the  90-deg  boundary,  so  the  reflection  coefficient  is  somewhat  higher 
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Figure  C3.  As  in  Figure  C2,  but  with  e  smeller  reflection  coefficient 
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Figure  C5.  A  case  with  modal  peaks  in  close  proximity,  broad  modal 
spreads,  low  noise,  and  moderate  reflection  coefficient 
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(X  —  0.43)  than  the  case  where  the  modes  are  well  isolated  (x  -  0.33).  As  can 
be  seen  in  Figure  C5 ,  variations  of  the  IKLE  estimate  about  the  true  spectrum 
do  not  have  a  large  effect  on  the  reflection  coefficient  because  both  modal 
regions  have  large  areas  and  the  noise  level  is  low.  Consequently,  the 
estimated  reflection  coefficient  is  only  in  error  by  about  2  percent. 

18.  If  the  noise  level  is  increased  for  the  case  shown  in  Figure  C5, 
error  estimates  remain  small,  partly  because  both  incident  and  reflected 
signals  remain  relatively  high.  Figure  C6  shows  a  case  where  N/S  —  0.1,  and 
all  other  conditions  are  the  same  as  in  Figure  C5.  The  IMLE  estimated 
spectrum  matches  the  true  spectrum  almost  perfectly,  and  the  error  in  the 
reflection  coefficient  is  less  than  1  percent.  Part  of  the  reason  for  this 
very  good  result  is  somewhat  fortuitous  in  that  the  structure  of  the  true 
spectrum  is  rather  well  matched  to  the  array  geometry.  Consequently,  fewer 
gages  are  lost  in  the  dynamic  elimination  process,  leaving  more  gages  with 
viable  information  to  resolve  the  true  spectral  shape.  In  this  case  11  gages 
are  used,  whereas  in  the  other  cases  shown  in  this  appendix,  only  6  to  8  gages 
are  used. 

19 .  If  the  amount  of  energy  in  the  reflected  mode  is  reduced  from  the 
case  shown  in  Figure  C6,  there  is  very  little  change  in  the  result.  Figure  C7 
shows  a  case  where  the  reflected  mode  has  only  2  percent  of  total  signal 
energy,  a  condition  where  the  reflection  coefficient  would  be  about  0.14  if 
the  modes  were  well  separated.  In  Figure  C7,  the  true  reflection  coefficient 
is  about  0.30,  mostly  because  the  tail  of  the  incident  mode  dominates  the 
calculation.  The  IMLE  estimate  remains  good,  with  the  error  in  reflection 
coefficient  being  less  than  3  percent. 

20.  Clearly,  Figures  Cl  to  C7  show  that  there  is  great  variety  in  the 
quality  of  reflection  coefficients  derivable  from  the  SAMSON  array  with  the 
IMLE  algorithm.  Results  depend  strongly  on  the  nature  of  the  overall  spectral 
shape,  the  amount  of  energy  in  the  reflected  mode,  and  the  amount  of  noise. 
There  is  also  a  likely  dependence  on  wave  frequency  because  short-wave  (high- 
frequency)  spectra  tend  to  be  better  resolved  than  long-wave  (low-frequency) 
spectra  with  an  array  of  finite  spatial  extent.  To  gain  a  broader  view  of  the 
ability  to  resolve  reflection  coefficients,  a  number  of  calculations  were 
done,  varying  the  test  frequency  and  all  six  of  the  parameters  governing  the 
shape  of  the  true  spectrum. 
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Figure  C6.  As  in  Figure  C5 ,  but  with  increased  noise 
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Figure  C7 .  As  in  Figure  CC,  but  with  less  energy  in  the  reflected  mode 
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21.  To  obtain  more  realistic  results,  ranges  of  values  for  background 
noise  were  approximated  from  the  SAMSON  spectral  parameters  listed  in  Appendix 
A.  As  discussed  in  paragraphs  38  to  41  in  the  body  of  this  report,  overall 
characteristic  wave  height  H,*,  was  found  from  the  total  volume  of  the 
frequency-direction  spectrum,  and  so  included  any  noise  that  was  present, 
whereas  the  wave  heights  that  characterized  incident  (H^,,  and  reflected 
(Hmo ,r)  energy  were  computed  after  first  removing  the  noise,  which  was  estimat¬ 
ed  by  averaging  the  smallest  10  values  of  the  integrated  direction  spectrum 
S(6a)  ,  defined  by  Equation  3  in  the  main  text.  Signal  variance  plus  noise  is 
then  16  times  the  square  of  ,  and  signal  variance  alone  is  estimated  as 

16  times  the  sum  of  the  squares  of  i  and  r  .  An  estimate  of  the 
noise-to-signal  ratio  is  thus  found  from 

N/S  .  (C12) 

16  <Ci+0 


or,  more  simply, 


N/S  = 


<o. 


Cr 


(C13) 


Figure  C8  shows  the  set  of  N/S  values  computed  from  Equation  C13  plotted  as 
a  function  of  .  For  moderate  and  high  ,  the  estimated  N/S  values 

in  Figure  C8  lie  mostly  in  the  range  0.01  to  0.1.  Significant,  but  not 
surprising,  is  that  N/S  tends  to  increase  for  small  wave  heights  (i.e.,  small 
signal).  For  H,^  <  0.5  m  ,  the  upper  limit  of  the  range  increases  to  about 
0.3.  While  the  numbers  shown  in  Figure  C8  are  only  approximate,  they  do 
suggest  the  order-of-magnitude  range  for  N/S  to  use  in  synthetic  data  tests. 

22.  Based  on  these  results,  tests  were  conducted  using  four  values  for 
N/S:  0.0001,  0.01,  0.1,  and  0.3,  with  the  first  one  being  a  very  low  value  to 
act  as  the  control  for  essentially  no  noise.  The  signal  part  of  the  synthetic 
spectrum  was  defined  by  varying  the  remaining  five  governing  parameters.  The 
parameter  set  consisted  of:  five  ?i  -  0.99,  0.98,  0.95,  0.90,  and  0.80;  four 
0p  i  -  0,  20,  40,  and  60  deg;  four  0p  r  -  120,  140,  160,  and  180  deg;  five 
L$t  -  10,  20,  30,  40,  and  50  deg;  and  four  A0r  -  10,  20,  30,  and  40  deg.  All 
combinations  of  all  parameters  were  employed,  resulting  in  1,600  different 
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Figure  C8.  Guidance  for  N/S  ratios  to  use  in  synthetic  spectra 

shapes  for  the  signal  spectra.  An  additional  dependency  on  wavelength  was 
also  examined  by  using  four  wave  frequencies  equal  to  0.06,  0.10,  0.14,  and 
0.21  Hz.  A  constant  water  depth  of  13  m  was  used  for  all  cases.  Altogether, 
25,600  test  cases  were  evaluated  for  error  in  reflection  coefficient. 

23.  Resulting  estimates  of  error  were  grouped  both  by  frequency  and  N/S 
ratio.  Groups  were  then  plotted  (Figures  C9  to  C12)  as  scatter  diagrams  of 
error  as  a  function  of  true  reflection  coefficient.  Results  for  0.06>Hz  waves 
are  shown  in  Figure  C9;  those  for  0.10-Hz  waves  are  shown  in  Figure  CIO;  those 
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for  0.14-Hz  waves  are  shorn  in  Figure  Cll;  and  results  for  0.21 -Hz  waves  are 
shown  in  Figure  C12.  A  few  (<  10)  extreme  outlier  errors  (near  400  percent) 
for  some  values  of  x  are  not  plotted  in  Figures  C9  to  C12  so  that  the  major 
point  groupings  could  all  be  plotted  at  the  same  scale.  Collectively,  Figures 
C9  to  C12  show  tendencies  for  the  error  in  reflection  coefficient  to  increase 
with  increasing  N/S,  increase  with  decreasing  wave  frequency,  and  increase 
with  decreasing  true  reflection  coefficient.  The  plots  also  show  a  strong 
tendency  for  the  algorithm  used  here  to  overpredict  the  true  reflection 
coefficient,  probably  because  the  noise  level  tends  to  be  underestimated,  as 
mentioned  in  the  discussion  of  Figures  C 3  and  C4,  above. 

24.  A  summary  of  the  results  of  this  analysis  is  obtained  by  computing 
means  and  standard  deviations  for  subsets  of  the  error  analysis  results 
(including  the  outlier  points  not  shown  in  Figures  C9  to  C12).  A  strong 
banding  of  data  is  seen  in  Figures  C9  to  C12  for  the  five  values  of  x 
corresponding  to  the  result  of  Equation  C5  for  the  five  PA  used  in  these 
tests.  Each  band  contains  about  250  data  points,  so  reasonably  representative 
means  and  standard  deviations  can  be  computed  for  these  five  x  •  Errors  for 
intermediate  values  of  x  follow  the  same  distribution  patterns  as  the  five 
main  clusters  in  all  of  the  figures,  so  the  summary  results  are  properly 
representative  of  the  whole  data  set. 

25.  Figure  C13  illustrates  the  summary  statistics.  Nominal  ranges  of 
frequency  and  noise- to-signal  ratio  in  which  reflection  coefficients  computed 
from  the  SAMSON  data  are  reasonably  reliable  can  be  deduced  from  this  figure. 
Both  the  bias  and  range  of  error  are  largest  for  small  x  at  all  noise 
levels .  The  high  range  of  error  suggests  that  results  are  highly  dependent  on 
specific  structures  of  the  true  spectra.  This  result  means  that  evaluating 
errors  for  low-x  SAMSON  data  would  have  to  be  done  on  a  case-by-case  basis 
because  the  spectral  shapes  vary  considerably.  Such  an  analysis  is  beyond  the 
scope  of  this  report,  but  it  seems  that  reflection  coefficients  near  0.1  or 
less  are  likely  to  be  considerably  in  error.  If  background  noise  levels  are 
high,  the  cutoff  for  reliable  reflection  coefficients  increases.  At 

N/S  -  0.3  ,  the  expected  (i.e.,  the  mean)  error  does  not  drop  below  20  percent 
until  x  is  about  0.3  or  higher. 

26.  For  low  to  moderate  noise  levels,  for  all  frequencies  tested,  and 
for  x  £  0.2  ,  the  expected  error  is  less  than  about  20  percent  and  the 
scatter  of  error  estimates  is  small,  suggesting  that  results  are  not  too 
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Reflection  Coefficient  Error  Statistics 
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dependent  on  specific  spectral  shape.  In  the  SAMSON  data,  estimated  noise 
levels  are  reasonably  low  for  greater  than  about  0.5  m.  Consequently, 

reported  reflection  coefficients  are  expected  to  be  good  at  the  20-percent 
level  for  x  >  0.2  and  H,*,  >  0.5  m.  An  uncertainty  of  a  factor  of  2  or  more 
can  be  expected  for  x  *  0.1  . 
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Appendix  D:  Notation 


dd 

df 

dff 

fP 

Vi 

fP,r 

H„o 

Hws.i 

^mo,r 


i 

k 


Coincident  spectral  density  between  gages  at  coordinates 
Xp,yp  and  xq,yq 

Mnemonic  for  two-digit  day  of  the  month 
Increment  of  frequency 
Increment  of  direction 

Element  of  directional  distribution  function  at  nth 
frequency  and  mth  direction 

nth  discrete  frequency  of  a  set  of  N  frequencies 
Spectral  peak  frequency 
Peak  frequency  of  incident  spectrum 
Peak  frequency  of  reflected  spectrum 
Spectrum-based  characteristic  wave  height 
Characteristic  height  of  incident  spectrum 
Characteristic  height  of  reflected  spectrum 


Complex  notation, 
Radian  wave  number 


m  Index  integer 

mm  Mnemonic  for  two-digit  representation  of  month  in  a 

calender  year 

M  Upper  limit  of  indices  denoted  m 

N  Upper  limit  of  indices  denoted  n 

N/S  Noise-to-signal  ratio 

Pt  Fraction  of  total  signal  variance  in  the  incident  mode 

of  a  bimodal  test  spectrum 

Qpq  Quadrature  spectral  density  between  gages  at  coordinates 

Xj,,yp  and  xq,yq 

RMS  Root -mean-square  difference  of  true  and  estimated 

spectra  expressed  as  a  percentage  of  SMJ< 
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RMS! 

RMSm 

s(f„) 

s(ea) 

S(fp.i.^-) 

S.(f«> 

Si 

Si(f„) 

Si<0 

s«* 

SM<«.) 

Sn 

S»(fp.i) 

Sp(fp,r) 

Sr 

sr(f„) 


Root -mean- square  difference  of  IMLE  spectrum  and  true 
spectrum  expressed  as  a  percentage  of 

Root -mean- square  difference  of  MLE  spectrum  and  true 
spectrum  expressed  as  a  percentage  of 

Frequency  spectral  density  at  nth  discrete  frequency 

Direction  spectral  density  at  mth  discrete  direction 

Frequency-direction  spectral  density  at  nth  discrete 
frequency  and  mth  discrete  direction 

Frequency -direction  spectral  density  at  peak  frequency 
of  the  incident  spectrum  and  at  mth  discrete  direction 

Frequency -direction  spectral  density  at  peak  frequency 
of  the  reflected  spectrum  and  at  mth  discrete  direction 

Auxiliary  frequency  spectral  density  at  nth  discrete 
frequency 

Estimate  of  total  incident  variance 

Incident  frequency  spectral  density  at  nth  discrete 
frequency 

Spectrum  at  mth  discrete  direction  estimated  with  the 
IMLE  algorithm 

Maximum  value  of  true  spectrum 

Spectrum  at  mth  discrete  direction  estimated  with  the 
MLE  algorithm 

Total  variance  in  noise  spectrum 

Estimated  total  variance  in  noise  spectrum  at  the  peak 
frequency  of  the  incident  spectrum 

Estimated  total  variance  in  noise  spectrum  at  the  peak 
frequency  of  the  reflected  spectrum 

Estimate  of  total  reflected  variance 

Reflected  frequency  spectral  density  at  nth  discrete 
frequency 

Spectral  peak  period 
Peak  period  of  incident  spectrum 
Peak  period  of  reflected  spectrum 
Mnemonic  for  two-digit  representation  of  year 
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yymmdd  Date  mnemonic,  a  concatenation  of  year,  month,  and  day 

mnemonics 

6$,  Directional  half-width  of  the  incident  signal  mode  of  a 

bimodal  test  spectrum 

69t  Directional  half-width  of  the  reflected  signal  mode  of 

a  bimodal  test  spectrum 

&8 1  Directional  width  parameter  of  the  incident  signal  mode 

of  a  bimodal  test  spectrum 

A 8r  Directional  width  parameter  of  the  reflected  signal  mode 

of  a  bimodal  test  spectrum 

Ax  Error  in  estimated  reflection  coefficient  as  a 

percentage  of  the  true  reflection  coefficient 

Axi  Error  in  reflection  coefficient  estimated  with  the  IMLE 

algorithm  as  a  percentage  of  the  true  reflection 

coefficient 

AXm  Error  in  reflection  coefficient  estimated  with  the  MLE 

algorithm  as  a  percentage  of  the  true  reflection 

coefficient 

A  Wavelength 

6m  Element  m  of  a  discrete  set  of  wave  directions 

6 p  Peak  wave  direction 

6p  y  Peak  direction  of  incident  wave  spectrum 

6 p  r  Peak  direction  of  reflected  wave  spectrum 

X  Reflection  coefficient  based  on  wave  heights 

Xj  Reflection  coefficient  of  a  spectrum  estimated  with  the 

IMLE  algorithm 

Xao  Reflection  coefficient  based  on  characteristic  wave 

heights 


XM 

Reflection  coefficient  of  a 
MLE  algorithm 

Xp.i 

Reflection  coefficient  at 
incident  wave  spectrum 

Xp,r 

Reflection  coefficient  at 
reflected  wave  spectrum 

spectrum 

estimated  with 

the 

• 

the 

peak 

frequency  of 

the 

the 

peak 

frequency  of 

the 

• 
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